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ABSTRACT 


V 


^  An  experimental  studyofj>oth  the  weakly  non-linear  as  well  as  the  three-dimensional 
nature  ofeboundary  layer  transition  is  conducted  using  the  active  surface  heating  tech¬ 
nique,  of  Liepmann  et  al.  In  the  present  study,' > this  technique  is  extended  to  provide  a 
means  for  controllably  and  repeatably  introducing  three-dimensional  disturbances  into 
a  laminar  boundary  layer;  'A  review  of  the  surface  heating  technique  is  presented  along 
with  a  discussion^ some  peculiarities  encountered  in  extending  this  technique  to  three- 
dimensiondT  geometries.  A  thorough  description  of  the  design  and  operation  of  a  pro¬ 
grammable  32-element  heater  array  and  the  supporting  instrumentation  are  also  given. 

/  The  heater  array  is  first  used  to  study  the  effect  of  weak  non-linearity  on  boundary 
/layer  transition.  By  keeping  the  forced  disturbances  as  two-dimensional  as  possible,  it  is 
/  shown  that  the  effects  of  weak  non-linearity  are  relatively  benign.  The  growth  rates  are 
seen  to  follow  the  linear  theory  up  to  perturbation  amplitudes  (r'w/r')  of  nearly  twelve 
percent.  The  only  deviation  from  the  linear  theory  arises  in  the  form  of  non-dinearly 
generated  harmonics  phase-locked  to  the  fundamental.  It  is  concluded  that  although 
these  non-Unearly  generated  harmonics  do  alter  the  wave  behavior  to  some  extent,  they 
are  by  themselves^ot  sufficient  to  explain  the  transition  from  small  linear  oscillations  to 
the  large  amplitude,  broadband,  three-dimensional  oscillations  characteristic  of  a  fully 
turbulent  boundary  layer. 


The  effect  of  threerdimensionality  on  boundary  layer  transition  is  then  investigated 
through  an  analytical  and  experimental  study  of  single  oblique  instability  waves.  This 
subject  has  remained  largely  unexplored,  as  such  disturbances  were  generally  thought 
to  be  more  stable  and  therefore  less  dangerous  than  their  two-dimensional  counterparts. 

•  /Through  a  series  of  experiments,  however,  it  is  shown  that  certain  conditions  exist  for 
which  oblique  waves  are  observed  to  be  more  unstable  than  any  two-dimensional  wave. 
It  is  shown  that  oblique  waves  exhibit  a  non-stationary  period-doubling  behavior  that 
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CHAPTER  1 


INTRODUCTION 


1.1  Historical  Background  and  Motivation 


The  subject  of  boundary  layer  transition  has  been  of  great  interest  for  over  fifty 


years  and  remains  a  very  active  area  of  research.  Historically,  the  most  fruitful  method 


by  which  transition  has  been  studied  is  a  consideration  of  the  stability  of  small  amplitude 


disturbances  of  the  form 


u(f,t)  =  u(y)e‘<“x-“‘> 


where  a  is  the  streamwise  wavenumber,  and  w  is  the  frequency  of  the  disturbance.  Both 


a  and  oj  are  in  general  complex.  Disturbances  of  the  form  1.1  are  most  commonly  referred 


to  as  Tollmien-Schlichting  waves  after  a  series  of  fundamental  papers  by  Tollmien  and 


Schlichting  from  1929  to  1940.  When  a  velocity  field  of  the  form  1.1  is  substituted 


into  the  linearized  equations  of  motion,  a  fourth  order  eigenvalue  problem  results  and 


is  typically  written  in  terms  of  the  vertical  component  of  the  disturbance  velocity  field 


v(y)  as 


-  2aV'  +  a4v  =  i*Re[(V  -  -)(t>"  -  a2v)  -  U"v 


U  =  U(y)  is  the  mean  velocity  profile,  and  Re  is  the  Reynolds  number.  Equation  1.2 


is  the  Orr-Sommerfeld  equation  derived  independently  by  Orr  (1907)  and  Sommerfeld 


(1908).  Solutions  are  easily  obtained  by  numerical  methods  and  are  usually  represented 


in  terms  of  the  neutral  stability  curve,  an  example  of  which  is  given  in  figure  1.1.  The 


existence  of  such  travelling  wave  disturbances  was  firmly  established  in  the  classical  ex¬ 


periments  of  Schubauer  and  Scramstad  (1947)  in  which  two-dimensional  oscillations  were 


artificially  introduced  into  a  laminar  boundary  layer  by  means  of  an  electromagnetically 


excited  vibrating  metal  ribbon.  Aided  by  the  extremely  low  turbulence  level  (0.02  Na¬ 


tional  Bureau  of  Standards’  wind  tunnel,  Schubauer  and  Scramstad  were  able  to  verify 


the  linear  theory  in  virtually  every  respect. 
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The  existence  of  the  disturbances  predicted  by  the  linear  theory  has  thus  been  estab¬ 


lished  beyond  doubt,  yet  the  connection  between  these  small  amplitude,  low  frequency, 


slowly  growing  oscillations  and  the  complicated,  three-dimensional  phenomenon  of  tur¬ 


bulence  still  remains  unclear.  Two  primary  directions  have  been  pursued  in  an  attempt 


to  relate  these  two  phenomena.  One  approach  is  a  consideration  of  non-linearity  in  the 


development  of  finite  amplitude  oscillations.  This  approach  has  been  taken  by  numerous 


investigators,  and  is  admirably  reviewed  in  a  recent  survey  by  Herbert  (1984).  Though 


the  results  of  the  various  analyses  vary  in  specific  details,  the  general  conclusion  reached 


by  this  line  of  investigation  is  that  non-linearity  alters  the  growth  rates  and  the  region  of 


unstable  parameter  space  only  slightly  and  cannot  in  and  of  itself  explain  the  transition 


from  linear  oscillations  to  fully  developed  turbulent  Row. 


Essential  to  a  more  complete  understanding  of  the  transition  process  is  some  account 


of  the  origin  of  three-dimensionality  in  the  flow.  As  early  as  1933,  Squire  related  the 


stability  of  three-dimensional  disturbances  in  the  form  of  oblique  waves  given  by 


*1 


to  the  well  understood  two-dimensional  problem.  By  a  simple  scaling  argument,  Squire 


showed  that  the  stability  of  oblique  waves  was  governed  by  an  equation  of  the  form 


1.2,  and  thus  the  eigenvalues  and  eigenvectors  could  be  obtained  for  this  type  of  three¬ 


dimensional  disturbance  just  as  in  the  two-dimensional  case.  The  principal  result  to 


come  from  this  analysis  was  the  observation  that  the  minimum  critical  Reynolds  number 


for  disturbances  of  the  form  1.3  occurs  for  the  purely  two-dimensional  case,  0  =  0.  It  is 


at  least  partly  for  this  reason  that  in  the  50  years  since  Squire’s  analysis,  very  few  experi¬ 


ments  exploring  the  stability  characteristics  of  oblique  waves  have  been  conducted.  Such 


disturbances  were  generally  considered  to  be  more  stable  and  therefore  less  dangerous 


than  two-dimensional  disturbances. 


It  was  not  until  the  experiments  of  Klebanolf,  Tidstrom,  and  Sargent  (1962)  that 
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interest  in  the  three-dimensional  nature  of  boundary  layer  transition  began  to  grow.  In 
those  experiments,  it  was  noticed  that  the  boundary  layer  tended  to  develop  a  spanwise 
modulation  in  the  disturbance  amplitude  with  a  seemingly  preferred  spanwise  wave¬ 
length,  A,.  This  behavior  was  enhanced  by  placing  physical  disturbances  on  the  surface 
of  a  flat  plate  at  pre-selected  spanwise  intervals.  A  very  regular  “peak-valley”  structure 
of  the  disturbance  waveform  was  then  clearly  seen.  The  “peaks”  correspond  to  regions 
of  maximum  wave  amplitude  and  the  “valleys”  to  regions  of  minimum  amplitude. 

Recent  flow  visualisation  studies  by  Saric  et  al.  (1984)  and  Kachanov  and  Levchenko 
(1984)  have  suggested  that  this  peak- valley  corrugation  of  the  boundary  layer  can  oc¬ 
cur  in  two  fundamentally  different  geometrical  patterns.  The  first  pattern,  which  was 
observed  by  Klebanoff  and  co-workers,  is  characterised  by  an  alignment  of  successive 
peaks  in  the  streamwise  direction;  i.e.  peaks  follow  peaks  and  valleys  follow  valleys. 
An  alternative  pattern  is  one  in  which  a  staggered  structure  develops;  i.e.  peaks  follow 
valleys  and  vise  versa. 

Numerous  mechanisms  have  been  proposed  in  an  attempt  to  explain  and  predict 
these  various  three-dimensional  wave  patterns.  Most  mechanisms  involve  a  weakly  non¬ 
linear  resonance  phenomenon  between  a  two-dimensional  Tollmien-Schlichting  wave  and 
some  form  of  three-dimensional  wave.  The  following  table  lists  several  of  the  resonance 
mechanisms  which  have  been  proposed: 
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Table  1.1  Various  Proposed  Resonance  Mechanisms 


Author 

Model 

Benny  and  Lin  (1960) 

Benny  (1961) 

Stuart  (1962,1971) 

-1-  B(cos/9z)e’(ax-"*‘) 

Herbert  and  Morkovin  (1980) 
Itoh  (1980) 

same  as  above,  but  to  3rd  order  in  A,B 

Craik  (1971) 

Raetz  (1959,1970) 

y|e«(ais-wtt)  +  get(aix+0a-v3t)  +  C’e»(a»*+/J«-w3t) 

c*l  —  c*j  +  Qrj  ,  u>i  =  a/j  +  wj 

Nayfeh  and  Bozatli  (1979) 

Agi(as-wt)  +  geHiax~tut)  _j_  Qgi(ax+0M-tut) 

Herbert  and  Morkovin  (1980) 

Ae'(«*‘“‘)  +  Beifia 

Benny  and  Gustavson  (1981) 

Aei(a,+fi.-»t)+  -squire  Modes” 

Each  model  differs  in  the  exact  form  of  the  three-dimensional  wave  chosen,  but  all 
predict  growth  rates  for  the  assumed  disturbances  that  are  in  excess  of  those  predicted 
by  the  linear  theory  alone.  The  so  called  "Squire  modes”  of  the  Benny-Gustavson  model 
are  solutions  of  the  linearized  equations  of  motion  of  the  form  1.3  but  with  zero  vertical 
component  of  the  disturbance  velocity  field,  i.e.  v{y)  =  0.  Such  disturbances  consist  of 
streamwise  and  spanwise  periodic  arrays  of  vertical  vortices  and  were  found  by  Squire 
(1933)  to  be  everywhere  stable.  They  are  variously  termed  “Squire  modes”  or  “vertical 
vorticity  modes”  in  the  literature. 


In  addition  to  the  above  mentioned  resonance  mechanisms,  several  analyses  have 
been  put  forth  treating  the  onset  of  three-dimensionality  as  a  secondary  instability  prob¬ 
lem.  Currently  receiving  much  attention  is  the  analysis  of  Herbert  (1983,1984)  which 
employs  Floquet  theory  to  study  the  instability  of  a  streamwise  periodic  base  flow  con¬ 
sisting  of  a  mean  flow  profile  plus  a  Tollmien-Schlichting  wave  that  is  neither  growing  nor 
being  damped.  This  is  strictly  valid  only  for  Pouiseiulle  flow,  but  is  suggested  to  provide 
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a  plausible  mechanism  for  the  boundary  layer  as  well.  This  analysis  is  appealing  in  that  it 
provides  theoretical  support  for  both  the  corrugated  three-dimensional  breakdown  seen 
by  Klebanoff  et  al.  as  well  as  the  staggered  structure  of  Saric  and  Kachanov. 

1.2  Objective  of  present  study 

Apart  from  those  studies  already  mentioned,  experimental  investigations  into  the 
three-dimensional  nature  of  boundary  layer  transition  remain  few  and  far  between.  At 
least  partly  responsible  for  the  lack  of  experimental  studies  is  the  fact  that  there  has  not 
existed  a  good  technique  for  introducing  three-dimensional  perturbations  into  a  laminar 
flow  in  a  controllable  and  repeatable  manner.  It  is  a  primary  goal  of  this  research  to 
describe  such  a  technique. 

The  technique  used  throughout  the  experiments  is  the  active  surface  heating  tech¬ 
nique  introduced  by  Liepmann,  Brown,  and  Nosenchuck  (1982).  A  review  of  the  surface 
heating  technique  will  be  presented  in  the  following  section  along  with  a  discussion 
of  some  peculiarities  introduced  by  the  use  of  such  a  technique  for  forcing  in  three- 
dimensional  geometries.  Chapter  2  will  then  describe  the  design  and  operation  of  the 
fully  programmable  32-element  heater  array  and  supporting  instrumentation  used  in  the 
present  experiments. 

A  second  goal  of  this  study  is  to  attempt  to  separate  as  much  as  possible  the  role  of 
non-linearity  from  the  role  of  three-dimensionality  in  the  rapid  breakdown  of  a  laminar 
boundary  layer  flow.  To  that  end,  chapter  3  will  be  devoted  entirely  to  a  study  of  two- 
dimensional,  non-linearly  developing  disturbances.  It  will  be  shown  that  non-linearity 
manifests  itself  primarily  in  the  generation  of  higher  harmonics  which  are  related  to  the 
self-interaction  of  the  fundamental  wave  with  itself.  It  will  be  further  shown  that  rather 
large  amplitude  two-dimensional  waves  can  be  easily  generated  when  three-dimensional 


effects  are  kept  to  a  minimum.  This  chapter  will  serve  the  dual  role  of  demonstrating 
the  phased  heater  array  capabilities  on  a  relatively  well  known  flow  (  2-D  T.S.  waves  ) 
while  at  the  same  time  establishing  a  reference  for  further  work  on  three-dimensional 


waves. 

As  three-dimensional  disturbances  are  considerably  more  complicated  and  have  been 
theoretically  and  experimentally  studied  in  much  less  detail,  we  restrict  our  attention  for 
the  present  study  to  one  particular  form  of  three-dimensional  disturbance;  that  of  oblique 
waves.  These  disturbances  are  perhaps  the  most  fundamental  type  of  three-dimensional 
disturbance,  as  they  are  normal  modes  of  the  linearized  equations  of  motion.  It  is  curious 
to  note  that  although  oblique  waves  are  widely  used  in  theoretical  studies  (as  can  be 
seen  in  table  1.1),  virtually  no  experimental  data  exists  on  this  type  of  disturbance.  In 
order  to  establish  a  framework  for  the  experimental  work,  Chapter  4  will  be  devoted  to 
a  detailed  discussion  of  some  of  the  unique  properties  of  oblique  waves.  It  will  be  seen 
that  while  many  similarities  between  oblique  waves  and  two-dimensional  waves  exist,  a 
number  of  stiking  differences  are  to  be  found  as  well. 

Finally,  Chapter  5  will  present  some  experimental  results  on  the  behavior  of  forced 
oblique  waves.  It  will  be  shown  by  careful  analysis  of  the  fluctuating  wall-shear  time 
series  and  of  the  disturbance  eigenfunctions  that  oblique  waves  exhibit  inherently  non¬ 
stationary  behavior.  This  behavior  will  be  analyzed  in  light  of  the  geometrical  structure 
revealed  in  Chapter  4,  and  will  be  shown  to  be  consistent  with  a  simple  mechanism 
which  has  been  known  to  play  a  fundamental  role  in  shear  flow  development  for  many 
years,  that  of  vortex  pairing.  The  relation  of  the  present  results  to  ideas  prevalent  in  the 
current  literature  will  be  discussed,  and  extensions  to  other  related  shear  flows  (wakes, 
shear  layers,  jets,  etc.)  will  be  made. 


1.3  Active  Surface  Heating  Technique 


Throughout  the  present  experiments,  the  active  surface  heating  technique  intro¬ 
duced  by  Liepmann  and  Nosenchuck  (1981)  was  used  to  excite  various  three-dimensional 
normal  modes  of  the  boundary  layer.  The  technique  is  based  on  the  viscosity-temperature 
relation  of  the  fluid.  Using  an  analysis  similar  to  that  of  Nosenchuck  (1982),  we  con¬ 
sider  the  zero  pressure  gradient,  z-momentum,  boundary  layer  equation  with  variable 
viscosity,  n  -  n{T)  -  n(y)  : 

du  du  du  1  dr  d2u  1  dfi  du  .  . 

-37 +  +«— = -— =  1/— -r —  (1.4) 

dt  dx  dy  p  dy  dy2  p  dy  dy 


Rearranging, 


du  du  .  .du 

+  u—  +  (v+Vtff)  — 


dt 


dx 


where  we  have  defined  an  “effective”  vertical  velocity, 


(1.5) 


v«//(y)  = 


ldpdT 
pdT  dy 


(1.6) 


As  it  is  well  known  that  the  boundary  layer  is  most  suceptable  to  disturbances  in  the 
vicinity  of  the  critical  layer  (see  for  example  Lin  (1955)  and  Meier  and  Maier  (1984)  ), 
we  would  like  to  evaluate  expression  1.6  at  the  critical  layer,  ycr.  We  therefore  need  an 
expression  for  the  temperature  profile  T[y)  above  the  heater. 


An  approximate  expression  for  the  temperature  profile  under  the  conditions  of 
steady,  parallel,  two-dimensional  boundary  layer  flow  can  be  derived  by  a  similarity 
analysis.  The  full  analysis  is  given  in  appendix  A,  but  for  the  present  only  the  result  is 
needed.  We  take  as  the  similarity  variable 


n  = 


y_ 

6p 


(1.7) 


where  (  is  the  streamwise  distance  measured  from  the  leading  edge  of  the  heater.  This 
is  shown  schematically  in  figure  A.l.  Above  the  heater,  the  temperature  profile  is  then 


given  by 


T(»j)  =  7„+0.776A7  ["  e-n*'*dti  (1.8) 

Jo 

where  Tw  is  the  temperature  at  the  wall,  and  AT  =  (7*,  -  Tw)  is  the  temperature 
difference  between  wall  and  freestream.  Combining  1.8  and  1.6,  we  obtain 

=  (...) 

It  is  desirable  to  express  the  temperature  difference  AT  in  terms  of  quantities  easily 
measurable  in  the  laboratory.  We  can  relate  AT  to  the  total  heat  flux  Q  f  introduced 
into  the  flow  by 

Qf  =  t>[  (1.10) 

Jo 

where  b  is  the  heater  span  and  L  the  streamwise  length  of  the  heater.  Evaluating 
qw  =  -kdT/dy  from  1.8,  equation  1.10  becomes 

Qf  =  O.mkATb(PrRclCf) 1/3  (111) 

It  may  be  noted  that  this  is  exactly  the  same  relation  that  would  be  obtained  through 
the  use  of  Lighthill’s  heat  transfer  formula. 

Qf  may  be  further  related  to  the  total  power  delivered  to  the  heater,  Qt  =  V2/  R 
Following  the  analysis  of  Taylor  (1986),  we  take 

Q»  =  <?/(1+^)  (1*2) 

where  Q,  is  the  heat  transferred  to  the  substrate  via  conduction.  This  conduction  loss 
can  be  estimated  by  assuming  that  the  length  scale  over  which  the  substrate  temperature 
gradients  are  maintained  will  scale  with  the  frequency  of  the  forcing,  w;  i.e.  the  relevant 
length  scale  in  the  substrate  is  9,  «s  (ic,/w)1/3  Thus  we  obtain 


Combining  with  1.11, 


^rkt(^TiPrRelCi  r'/s  (,u) 

Under  the  present  conditions  this  quantity  is  of  order  10~l,  so  we  may  approximate  1.12 
as  simply  Qf  «  .9 Qt.  Combining  this  result  with  1.9  and  1.11  and  simplifying,  we  obtain 
the  final  approximate  result  that 


V'fftitV) 


I  Q*  c-r„v3 /9wte 
pdT  kb(L*t)l,i 


(1.15) 


This  expression  gives  the  effective  vertical  velocity  introduced  by  the  heating  tech¬ 
nique.  It  is  in  a  convenient  form  for  pointing  out  several  features.  First,  it  is  noted  that 
the  fluid  response  depends  linearly  on  the  power  input,  Qt.  This  is  in  agreement  with 
experimental  results  as  can  be  seen  in  figure  3.4  in  which  the  disturbance  amplitude 
measured  far  downstream  of  the  heater  is  seen  to  be  a  very  nearly  linear  function  of  the 
heater  power.  In  this  case,  there  would  be  additional  multiplicative  factors  introduced 
by  the  integrated  disturbance  growth  and  the  conversion  to  wall  shear,  but  the  linear 
dependence  upon  Qt  remains  the  same.  This  is  very  convenient  to  know,  for  it  means 
that  in  order  to  increase  the  disturbance  amplitude  by  say  a  factor  of  2,  one  simply  needs 
to  increase  the  heater  power  by  a  factor  of  2. 


Another  feature  of  equation  1.15  is  that  the  vertical  dependence  of  the  effective 
velocity  is  now  made  explicit.  For  instance,  evaluating  1.15  at  the  wall  and  integrating 
over  (,  we  obtain 


w«//l 


»=o 


1  dfi  Qt 
p  dT  kbL 


(1.16) 


This  is  easily  shown  to  be  identical  to  the  expression  for  the  effective  wall  velocity  as 
found  in  Nosenchuck.  Nosenchuck  then  argued  that  the  streamwise  extent  of  the  heater 
should  be  wide  compared  to  the  local  boundary  layer  thickness  to  allow  for  penetration 
of  the  heat  flux  to  the  critical  layer.  The  nature  of  this  penetration  is  seen  directly 


IP 
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in  1.15.  At  the  wall,  the  effective  velocity  perturbation  decreases  as  £-*/3  due  to  the 
increasing  thermal  boundary  layer  thickness  and  the  correspondingly  decreasing  wall 
temperature  gradient.  Away  from  the  wall,  however,  the  exponential  term  domimates, 
and  the  effective  velocity  now  increases  exponentially  with  £  as 

".//(£.»)  (1-17) 

As  we  move  further  away  from  the  wall,  this  spatial  localization  of  the  effective  velocity 
increases  greatly,  as  the  streamwise  dependence  on  £  is  multiplied  by  a  cubic  term  in  y. 
Therefore  in  the  vicinity  of  the  critical  layer,  the  effective  velocity  increases  exponentially 
with  streamwise  distance  along  the  heater  reaching  a  maximum  at  the  heater  trailing 
edge.  Though  our  similarity  analysis  extends  only  to  the  heater  trailing  edge,  it  is  easily 
seen  that  a  further  increase  in  £  will  result  in  a  decrease  in  the  effective  velocity  as  there 
is  no  longer  any  heat  input  at  the  wall,  and  all  existing  thermal  gradients  will  decay 
rapidly. 

In  the  present  study,  the  streamwise  extent  of  the  heater  is  7.6  mm  and  the  dis¬ 
placement  thickness  6’  is  typically  on  the  order  of  1  mm.  The  heater  is  therefore  rather 
large  compared  to  the  boundary  layer  thickness.  Also  of  interest  is  the  wavelength  of 
the  disturbance  which  in  the  present  experiments  is  on  the  order  of  25  mm.  The  heater 
is  seen  to  cover  a  significant  fraction  of  this  distance  as  well.  Due  to  the  exponential 
dependence  on  £  as  seen  in  equation  1.17,  however,  the  spatial  localization  of  the  forcing 
is  maintained,  and  the  resulting  resolution  of  the  forced  oscillations  in  both  time  and 
space  is  excellent. 


As  the  discussion  thus  far  has  dealt  with  a  two-dimensional  heater,  it  is  necessary 
to  make  a  few  comments  with  regard  to  the  spanwise  segmentation  of  the  present  heater 
array.  As  will  be  shown  in  the  next  chapter,  the  present  heater  array  consists  of  32 
elements  which  are  separated  by  a  heater-to- heater  gap  of  0.76  mm.  This  gap  is  small, 
but  it  is  not  negligible  when  compared  to  the  displacement  thickness  or  critical  layer 


i&S&SS 


height  both  of  which  are  of  the  order  of  1  mm.  The  input  disturbance  cannot  then  be 
considered  as  a  single  eigenmode,  but  must  rather  be  represented  as  a  sum  of  normal 
modes  over  all  possible  spanwise  wavenumbers  with  the  frequency  of  oscillation  set  by 
the  heater  forcing.  Thus  we  have  as  our  input  disturbance 

u(*,t)  =  (US) 

0 

Though  no  measurements  were  made  of  this  transient  disturbance  development,  it  is 
believed  that  the  preferential  amplification  of  the  boundary  layer  for  disturbances  of  low 
spanwise  wavenumber  will  insure  that  the  dominant  mode  amplified  will  indeed  be  the 
lowest  mode.  This  preferential  amplification  is  seen  in  figure  5.11  which  is  taken  from 
the  calculations  of  Mack  (1984).  Thus  the  high  wavenumber  contributions  to  the  sum  in 
equation  1.18  will  be  strongly  damped,  and  the  resulting  disturbance  waveform  observed 
far  downstream  will  be  very  close  to  a  single  normal  mode. 


1.4  On  Comparisons  with  Other  Experiments 

In  the  present  study,  the  primary  flow  diagnostics  are  hot-fiim  wall  shear  sensors,  and 
the  results  are  correspondingly  presented  in  terms  of  normalized  fluctuating  wall  shear, 
t'w/Tw  In  the  current  literature,  however,  hot-wire  anemometry  is  the  most  commonly 
used  method  of- measurement,  and  results  are  typically  presented  in  terms  of  velocity 
fluctuations  at  the  critical  layer  normalized  by  the  free  stream  velocity,  umoac /(/<».  It  is 
therefore  necessary  to  relate  these  two  quantities  in  order  that  reasonable  comparisons 
between  experiments  may  be  made. 

Close  to  the  wall,  a  reasonable  approximation  to  the  form  of  the  disturbance  eigen¬ 
function,  u(y),  can  be  taken  as 

«(y )  *  Uma*«n(«vy)  (1.19) 
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A s  the  maximum  in  the  disturbance  eigenfunction  occurs  very  nearly  at  the  critical  layer 
yer,  we  can  evaluate  the  vertical  wavenumber  kv  in  terms  of  ycr  as 


f  —  . 

9  2ye 


(1.20) 


Assuming  that  this  expression  for  the  vertical  wavenumber  holds  at  the  wall,  we  can 
evaluate  the  fluctuating  wall  shear  as 


,  du, 

r*  =^U  = 


dy'v=°  ~  2yer 


(1.21) 


Normalizing  tJ  by  the  mean  wall  shear  fm  and  umaz  by  the  freestream  velocity  (J^, 

*’•*  U max 


rm  2fwyer  Ua 


Using  the  Blasius  value  for  the  mean  wall  shear, 


(1.22) 


we  obtain 


rw  =  0.332 pU*Re,-l/i 


rw' /f„  _  KvRtxl/i 

Um**/f^oo  0.664f/ooycr 


(1.23) 


(1.24) 


Using  the  approximation  that  y er  »  6* ,  and  that  Rtf  =  1.721  Rtx1^,  we  obtain 


rj/fw 


m  2.75 


(1.25) 


U  max/^o 

This  is,  of  course,  only  a  rough  estimate,  but  it  provides  us  with  an  approximate  figure 
for  comparing  the  magnitudes  of  wall-shear  fluctuations  with  velocity  fluctuations. 

This  approximate  relation  can  be  verified  experimentally  by  simultaneously  measur¬ 
ing  these  two  quantities.  The  full  details  of  the  experimental  aparatus  will  be  discussed 
in  the  next  chapter.  For  the  present,  though,  it  is  necessary  only  to  note  the  result.  A 
cylindrical  hot-film  velocity  probe  was  positioned  at  the  maximum  in  the  eigenfunction 
profile.  In  this  case,  that  corresponds  to  a  distance  of  0.3  mm  above  the  plate  surface. 
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The  velocity  probe  was  positioned  directly  above  a  hot-film  wall-shear  probe  that  was 
flush  mounted  on  the  plate  surface.  A  two-dimensional  Toilmien  Schlichting  wave  was 
then  generated.  Figure  1.2  shows  the  simultaneously  recorded  time  traces.  The  max¬ 
imum  amplitude  of  the  fluctuating  velocity  umax/U ^  is  seen  to  be  approximately  two 
percent,  while  the  peak  amplitude  of  the  fluctuating  wall-shear  signal  rj/tv  is  approxi¬ 
mately  five  percent.  Thus,  the  experimentally  measured  ratio  of  2.5  compares  quite  well 
with  the  predicted  value  of  2.75. 


CHAPTER  2 


EXPERIMENTAL  FACILITY  AND  INSTRUMENTATION 
2.1  GALCIT  High  Speed  Water  Tunnel 

All  of  the  experiments  in  the  present  study  were  conducted  in  the  GALCIT  (Grad¬ 
uate  Aeronautical  Laboratories,  California  Institute  of  Technology)  High  Speed  Water 
TVinnel  (HSWT).  A  complete  description  of  the  HSWT  can  be  found  in  Ward  (1976). 
As  shown  in  figure  2.1,  the  HSWT  is  a  closed-circuit  facility.  The  test  section  was  of 
circular  cross  section  with  an  internal  diameter  of  35.6  cm  and  an  overall  streamwise 
length  of  1.05  m.  The  test  section  is  preceded  on  the  upstream  end  by  a  20:1  contraction 
section  that  is  three  meters  long,  and  it  is  followed  on  the  downstream  side  by  a  6:1 
diffusor  section  that  is  six  meters  in  length.  A  stainless  steel  honeycomb  and  a  40-mesh 
screen  are  positioned  at  the  beginning  of  the  contraction  section.  Throughout  the  present 
experiments,  the  tunnel  pressure,  which  is  adjustable  over  the  range  0.2  to  6  bar,  was 
maintained  near  atmospheric  pressure. 


A  35  cm  wide  lucite  flat  plate  was  mounted  in  the  test  section,  completely  spanning 
the  test  section  diameter.  Throughout  the  experiments,  the  freestream  velocity  was 
maintained  at  approximately  120  cm/s  resulting  in  a  Reynolds  number  of  1.4  x  10°  based 
upon  the  total  plate  length.  The  flow  quality  was  generally  found  to  be  quite  acceptable. 
Using  a  cylindrical  hot-film  probe,  the  r.m.s.  turbulence  intensity  was  measured  to  be 
0.05  percent  in  the  freestream. 


2.2  Flat  Plate  Model 

Figure  2.2  shows  the  flat  plate  test  model  used  throughout  the  experiments.  The 
plate  was  machined  from  lucite  and  consisted  of  a  central  flat  section  preceded  by  a 
6:1  elliptical  leading  edge  and  followed  by  an  adjustable  trailing  edge  flap.  The  overall 


plate  dimensions  (including  leading  and  trailing  edges)  were  104  cm  in  length  by  35 
cm  in  width  by  2.0  cm  in  thickness.  The  leading  edge  was  chosen  to  be  a  6:1  ellipse,  as 
that  geometry  provided  a  reasonable  laminar  flow  profile  over  a  large  portion  of  the  plate 
center  body  while  maintaining  a  small  but  finite  region  of  adverse  pressure  gradient  close 
to  the  leading  edge.  The  phased  heater  array  was  deliberately  located  in  this  region  so 
as  to  take  advantage  of  the  enhanced  growth  rates  offered  by  the  adverse  gradient  and 
thus  augment  the  effectiveness  of  the  heating  technique. 

The  central  flat  section  of  the  plate  contained  most  of  the  forcing  and  flow  diagnostic 
equipment.  The  phased  heater  array  was  flush  mounted  7.0  cm  from  the  leading  edge. 
At  a  typical  freestream  velocity  of  120  cm/s,  the  heater  Reynolds  number  Res-  was  ap¬ 
proximately  520.  Thus  with  respect  to  the  neutral  stability  curve,  the  heater  was  located 
very  close  to  the  minimum  critical  Reynolds  number.  When  not  in  operation,  the  heater 
presents  a  very  minimal  disturbance  to  the  flow.  For  example,  non-dimensionalizing  by 
the  wall  friction  velocity  u*  =  (rw/p)1^3,  the  protrusion  of  the  heater  into  the  flow  is 
seen  to  be  a  mere  y+  =  u't/u  »  1.5,  where  t  is  the  thickness  of  the  etched  copper  layer. 

Hot-film  wall-shear  sensors  were  located  at  various  streamwise  and  spanwise  loca¬ 
tions  far  downstream  of  the  heater  array,  and  provided  the  primary  flow  diagnostics. 
The  plate  also  contained  a  static  pressure  tap  which  together  with  a  pitot  probe  located 
upstream  of  the  contraction  section  was  used  to  obtain  the  freestream  flow  conditions 
and  thereby  calibrate  the  wall-shear  sensors.  For  several  experiments,  a  dye  injector  was 
also  located  on  the  plate  for  flow  visualization. 

All  electrical  leads  and  tubing  were  connected  to  the  bottom,  non-active  side  of  the 
plate  and  exited  through  sealed  plugs  on  the  bottom  wall  of  the  test  section.  As  the 
substantial  number  of  connections  together  with  the  finite  thickness  of  the  plate  itself 
provided  a  non-negligible  flow  blockage,  the  trailing  edge  flap  was  used  to  position  the 
leading  edge  stagnation  point  such  as  to  minimize  the  overall  plate  pressure  gradient. 


2.3  Flow  Documentation 


As  mentioned  previously,  the  plate  contained  a  6:1  elliptical  leading  edge.  This 
created  a  region  of  adverse  pressure  gradient  immediately  downstream  of  the  leading 
edge.  From  static  pressure  measurements  on  earlier  models  with  similar  leading  edge 
configurations,  it  was  found  that  the  flow  returned  to  an  essentially  undisturbed  Blasius 
flow  very  quickly.  In  order  to  measure  the  mean  velocity  profile,  a  cylindrical  hot- 
film  velocity  probe  (TSI  model  1260-10W)  was  used.  The  supporting  structure  for  the 
probe  was  designed  with  a  streamlined  cross-section  so  as  to  provide  adequate  rigidity 
while  minimizing  flow  induced  vibration  due  to  vortex  shedding.  The  boundary  layer 
thickness  699%  at  the  probe  location  was  approximately  3.0  mm,  and  the  probe  traversed 
this  distance  in  32  increments  of  0.1  mm  per  step.  The  probe  was  positioned  by  a  stepper 
motor  which  was  interfaced  to  the  central  microprocessor.  A  physical  stop  was  provided 
which  limited  the  probe  to  a  minimum  vertical  position  of  0.125  mm  (y+  ss  3.5). 

The  measured  boundary  layer  profile  is  shown  in  figure  2.3  with  the  computed 
Blasius  curve  included  for  comparison.  As  can  be  seen,  the  measured  profile  at  this 
x-location  is  very  nearly  Blasius.  Thus,  the  local  behavior  of  the  forced  disturbances 
as  measured  at  locations  far  from  the  leading  edge  should  agree  quite  well  with  the 
predictions  of  the  linear  theory  for  a  Blasius  mean  profile.  Measurements  of  integrated 
quantities  such  as  the  disturbance  amplitude,  however,  will  depend  upon  the  exact  pres¬ 
sure  gradient  history  encounterred  by  the  disturbance  as  it  travels  downstream  from  the 
heater  to  the  point  of  measurement.  For  this  reason,  most  results  in  this  study  will 
emphasize  measurements  of  local  quantities  as  opposed  to  integrated  quantities.  In  this 
way,  comparisons  with  other  experiments  or  numerical  calculations  may  be  made. 
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2.4  Phased  Heater  Array 

Figure  2.4(a)  shows  a  planform  view  of  the  heater  array.  The  array  consists  of 
32  individual  heater  elements,  each  with  one  lead  running  to  the  side  of  the  array  for 
connection  to  the  driving  electronics  and  the  other  lead  connected  to  a  common  ground. 
An  enlarged  view  of  a  single  heater  is  shown  in  figure  2.4(b).  The  overall  dimensions 
of  each  individual  element  are  1.0  cm  in  width  by  0.76  cm  in  the  streamwise  direction. 
Within  each  heater  the  line  width  is  0.254  mm,  and  the  total  line  length  is  approximately 
16.3  cm.  This  results  in  a  nominal  heater  resistance  of  0.50  ohms.  The  low  heater 
resistance  was  chosen  in  order  to  enable  large  power  input  into  the  flow  at  relatively 
low  voltages.  The  advantages  of  this  scheme  are  twofold.  First,  the  low  heater  voltages 
prevented  the  occurrence  of  electrolysis  and  the  associated  bubble  formation.  Secondly, 
and  of  equal  importance,  was  the  fact  that  very  large  amounts  of  power  could  be  easily 
supplied  at  such  low  voltages  by  simply  using  two  12  volt  batteries  as  the  sole  power 
supply.  Separating  each  heater  is  a  gap  of  0.76  mm  making  the  entire  array  34  9  cm 
in  width  which  spans  the  GALCIT  High  Speed  Water  Tunnel  cross  section  of  35.5  cm. 
The  array  is  fabricated  from  a  commercially  available  material  consisting  of  a  0.035  mm 
layer  of  copper  deposited  on  an  epoxy-glass  laminated  substrate  and  was  etched  using 
conventional  photolithographic  techniques.  All  design  and  construction  was  carried  out 
at  GALCIT. 


2.5  Instrumentation 


To  control  and  supply  power  to  the  heater  array,  a  good  deal  of  supporting  instru¬ 
mentation  was  necessary.  Figure  2.5  shows  a  schematic  diagram  of  the  control  appara¬ 
tus.  The  entire  array  was  under  control  of  an  8-bit  microprocessor  (Intel  8085)  which 
established  the  frequency,  relative  phase,  and  amplitude  of  the  32  heating  elements.  A 


Kp 


programmable  frequency  generator  was  employed  to  generate  sinusoidal  signals  covering 
the  range  from  0  to  50  Hz.  The  frequency  was  programmable  in  256  discrete  steps  with 
an  incremental  step  size  of  less  than  0.2  hz.  As  the  heat  flux  into  the  flow  Q /  is  pro¬ 
portional  to  the  square  of  the  applied  heater  voltage  (or  current),  the  available  forcing 
frequency  range  seen  by  the  flow  was  actually  twice  that  of  the  input  or  from  0  to  100 
Hz.  This  completely  spanned  the  range  of  unstable  Tollmein-Schlichting  frequencies  at 
the  flow  conditions  encountered  throughout  the  experiments. 

To  produce  three-dimensionality  in  the  form  of  arbitrary  spanwise  phase  distribu¬ 
tions,  a  bank  of  32  individually  programmable,  unity-gain  phase  shift  stages  was  pro¬ 
vided.  These  circuits  were  again  under  control  of  the  central  microprocessor  and  were 
programmable  in  16  discrete  increments.  The  minimum  heater  to  heater  phase  shift  was 
5.6*  and  the  maximum  possible  was  90*.  The  phase  control  was  primarily  used  for  the 
forcing  of  single  oblique  eigenmodes,  and  it  was  thus  found  convenient  to  tie  the  output 
of  phase  shift  n  -  1  to  the  input  of  phase  shift  n  resulting  in  a  linear  phase  shift  across 
the  span  of  the  array. 

For  example,  with  each  heater  programmed  for  an  11.2°  phase  shift,  a  total  phase 
shift  of  360*  would  result  over  the  35  cm  span  of  the  array.  Thus,  a  disturbance  would 
he  generated  with  a  spanwise  wavelength  A.  of  35  cm  Under  the  present  experimental 
conditions,  a  typical  streamwise  wavelength  Ax  would  be  on  the  order  of  2  cm,  and  the 
resulting  disturbance  would  therefore  be  an  oblique  eigenmode  with  wavenumber  vector 
ii  oriented  at  an  angle  9  =  tan-1(A*/A,)  «s  4°  from  the  freestream  flow  direction. 

In  addition  to  the  phase-shift  capability,  the  array  control  electronics  also  allowed 
for  amplitude  modulation  of  each  individual  heater  voltage.  As  with  the  phase  shift 
circuitry,  these  were  again  under  control  of  the  microprocessor  with  outputs  ranging 
from  0.5  volts  to  8.0  volts  in  increments  of  0.5  volts.  These  were  used  to  create  arbitrary 
amplitude  distributions  across  the  span  of  the  array.  As  an  example,  the  array  could  be 


t.  tt.  |l .  |r.  t».  jt, 


n.  ij.  a».  |i,  n, 


i—  |»  >■»  *  *_>  »  a  a  t  ■.*! 


-  19- 


programmed  for  a  sinusoidal  span  wise  amplitude  distribution  producing  a  disturbance 


waveform  very  similar  to  that  produced  by  the  vibrating  ribbon  technique. 


Finally,  before  reaching  the  heater  array,  each  signal  went  through  a  power  amplifi¬ 


cation  stage.  There  were  32  such  circuits,  one  to  supply  power  to  each  heating  element. 


The  circuit  as  seen  in  figure  2.6  was  a  simple  push-pull  emitter  follower  with  Darlington 


power  transistors  to  supply  the  necessary  current  amplification  {kps  «  4000).  The  push- 


pull  configuration  was  used  in  order  to  allow  for  symmetrical  forcing  about  ground  and 


thus  minimize  the  unnecessary  D.C.  heating.  During  the  experiments,  the  output  voltage 


was  limited  to  5.0  volts  resulting  in  a  maximum  r.m.s.  power  output  of  approximately 


25  watts/heater  or  over  700  watts  for  the  entire  array. 


In  order  to  obtain  a  physical  picture  of  the  types  of  disturbances  which  can  be 


produced  with  the  array,  three  separate  cases  are  considered.  Figure  2.7  shows  a  three¬ 


dimensional  perspective  view  of  the  resulting  Tollmein-Schlichting  wave  surfaces  for  each 


case.  In  each  plot  the  axes  are  time,  spanwise  distance,  and  the  normalized  wall-shear 


amplitude,  r'w/fu.  To  obtain  these  plots,  the  programmability  of  the  array  was  em¬ 


ployed.  Rather  than  move  the  sensor  relative  to  a  fixed  disturbance  as  is  customarily 


done,  the  disturbance  in  this  case  was  moved  relative  to  a  fixed  sensor  located  on  the 


plate  centerline  at  a  Reynold’s  number  of  1240.  In  this  way,  the  measurement  apparatus 


is  greatly  simplified. 


Figure  2.7(a)  shows  a  two-dimensional  Tollmien-Schlichting  wave  surface.  In  this 


case,  each  heater  was  forced  with  the  same  phase  and  amplitude  across  the  span  of 


the  array.  The  two-dimensionality  of  the  waves  has  been  measured  at  several  spanwise 


locations  and  is  seen  to  be  excellent.  Using  this  type  of  two-dimensional  forcing,  wave 


amplitudes  of  15  percent  have  been  obtained.  Figure  2.7(b)  shows  a  single  oblique 


eigenmode.  The  array  in  this  case  was  programmed  for  a  constant  phase  shift  from 


heater  to  heater  while  maintaining  a  constant  amplitude  across  the  span  of  the  array. 
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Though  such  oblique  waves  have  been  known  to  be  unstable  for  over  50  years,  they  have 


never  been  systematically  studied  due  to  lack  of  an  appropriate  experimental  technique. 


Finally,  figure  2.7(c)  shows  a  wave  packet;  i.e.  a  disturbance  localized  in  both 


space  and  time.  The  forcing  in  this  case  consisted  of  4  cycles  of  a  sinewave  forced  by  a 


single  heater  element  with  all  other  elements  turned  off.  The  single  active  element  was 


then  swept  across  the  array  creating  the  spanwise  disturbance  seen  in  the  figure.  These 


plots  are  shown  at  the  present  time  merely  to  give  a  qualitative  picture  of  the  possible 


types  of  disturbances  that  can  be  generated  by  the  phased  heater  array.  A  much  more 


detailed  look  at  the  experimentally  generated  wave  patterns  will  be  given  in  chapters  3 


(two-dimensional  waves)  and  5  (oblique  waves). 


2.6  Flow  Diagnostics  and  Data  Acquisition 


This  section  describes  the  flow  diagnostic  equipment  used  throughout  the  experi¬ 


ments.  The  primary  means  of  measurement  was  provided  by  hot-film  wall-shear  sensors 


(TSI  model  1240)  that  were  flush  mounted  on  the  plate  at  various  streamwise  and  span- 


wise  locations.  The  sensor  dimensions  were  1.0  mm  by  0.125  mm  (aspect  ratio  8:1)  and 


were  oriented  such  as  to  measure  the  streamwise  component  of  the  wall  shear.  These 


probes  were  used  to  obtain  measurements  of  growth  rates,  phase  velocities,  spanwise 


phase  relationships,  etc.  Measurements  could  be  made  at  streamwise  locations  separated 


by  as  little  as  3  mm,  a  distance  corresponding  to  roughly  0.1  Az,  where  Az  is  a  typical 


streamwise  Tollmien-Schlichting  wavelength.  In  this  way,  measurements  of  instantaneous 


growth  rates  could  be  obtained  even  when  the  disturbances  exhibited  non-linear  or  non¬ 


stationary  behavior.  It  is  important  to  emphasize  the  non-stationary  aspect  at  this  point. 


If  the  wave  growth  were  purely  stationary,  then  a  single  probe  could  be  traversed  in  the 


streamwise  direction  to  obtain  the  disturbance  growth  rates.  As  soon  as  the  problem 
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becomes  non-station  ary,  however,  measurements  at  different  locations  must  necessarily 
be  made  simultaneously  in  order  to  capture  the  relevant  phenomena. 

Both  the  wall-shear  sensors  and  the  velocity  sensor  were  operated  in  the  constant 
temperature  mode.  The  probe  resistance  (and  therefore  temperature)  was  maintained 
at  a  constant  value  by  the  bridge  circuit  of  figure  2.8.  The  circuit  is  a  simple  voltage 
feedback  loop  with  an  emitter-follower  (2N2270)  to  supply  the  necessary  current.  A  200 
ohm  potentiometer  was  used  to  adjust  the  probe  overheat,  A T/T  =  [R\0t  -  Rcoid) / Rcoid, 
which  was  typically  on  the  order  of  6  percent.  The  bridge  circuit  also  contained  a  current 
limiting  resistor  which  maintained  the  probe  current  at  a  safe  operating  level.  The  output 
from  the  top  of  the  bridge  was  then  buffered  and  low-pass  filtered  at  200  Hz  to  prevent 
aliasing  in  data  acquisition.  Power  was  supplied  by  two  12-volt  marine  batteries  which 
eliminated  the  60  Hz  line  noise  usually  introduced  by  rectifying  power  supplies.  The 
noise  level  was  thus  kept  to  a  quite  acceptable  low  level  of  5  millivolts  peak-to  peak. 
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CHAPTER  3 

TWO  DIMENSIONAL  WAVE  EXCITATION 
3.1  Introduction 

In  this  chapter,  we  will  examine  the  results  from  several  experiments  on  the  forcing 
of  two-dimensional  Tollmien-Schlichting  waves.  Emphasis  will  be  put  on  understanding 
the  nature  of  the  non-linear  effects  arising  as  the  disturbance  amplitude  is  increased 
As  indicated  earlier,  numerous  studies  have  been  undertaken  to  explore  the  non-linear 
aspects  of  boundary  layer  transition.  The  question  thus  arises  as  to  why  the  need  for 
another  study  on  this  same  subject.  The  reason  for  the  present  investigation  is  two-fold. 

First,  the  present  state  of  knowledge  on  the  subject  is  far  from  complete.  As  a 
specific  example,  consider  the  following  question:  at  what  disturbance  amplitude  does 
the  linear  theory  cease  to  hold?  That  is,  we  seek  the  threshold  amplitude  beyond  which 
either  non-linearity  or  three-dimensionality  or  both  must  be  taken  into  account  in  order 
to  explain  the  disturbance  behavior.  This  specific  point  provides  a  convenient  bench¬ 
mark  by  which  experiments  and  theory  may  be  compared,  as  the  threshold  amplitude 
is  a  frequently  quoted  figure  in  the  current  literature.  Furthermore,  from  a  practical 
standpoint,  the  threshold  amplitude  is  significant  in  that  it  is  a  measure  of  the  point 
at  which  large  dynamical  growth  rates  take  over  from  the  much  smaller  viscous  growth 
rates  predicted  by  the  linear  theory.  It  therefore  gives  an  approximate  measure  of  the 
point  at  which  transition  to  turbulence  may  be  said  to  occur. 

A  significant  number  of  investigators  report  a  threshold  amplitude  on  the  order 
of  one  percent;  i.e.  at  u'ma*/£A»  **  1%,  a  strong  departure  from  the  linear  theory  is 
observed.  This  figure  is  reported  in  the  experiments  of  Saric,  Kozlov,  and  Levchenko 
(1984)  and  Kachanov  and  Levchenko  (1984),  the  theoretical  analyses  of  Craik  (1971) 
and  Herbert  (1983),  and  the  numerical  simulations  of  Spalart  (1984). 


In  contrast  to  these  investigations  are  several  experimental  studies  which  report 
no  apparent  threshold  amplitude  for  disturbances  of  5  to  6  percent.  Strykowski  and 
Sreenivasan  (1985),  for  example,  using  a  purely  two  dimensional  mechanical  suction  and 
blowing  device  report  wave  amplitudes  of  close  to  6  percent.  The  agreement  between 
experimental  data  and  the  linear  eigenfunction  are  good.  Williams,  Fasel,  and  Hama 
(1984)  likewise  report  waves  of  5.5  percent,  again  with  excellent  agreement  to  the  linear 
theory.  The  forcing  mechanism  in  this  case  was  a  small  wire  located  at  the  critical  layer 
oscillated  by  a  mechanical  lever-cam  arrangement. 

-  Clearly  there  is  some  disagreement  on  this  issue.  It  is  to  be  noted  that  in  all 
of  the  studies  reporting  a  threshold  amplitude  near  one  percent,  the  departure  from 
linearity  is  coincident  with  the  onset  of  three-dimensionality.  This  is  evident  in  the  flow 
visualization  of  Saric  et  al.  and  Kachanov  et  al.  as  well  as  the  numerical  particle  tracing 
of  Spalart.  On  the  other  hand,  the  experiments  of  Williams  et  al.  clearly  document  the 
fact  that  the  observed  waves  in  their  experiment  are  two-dimensional.  This  discrepancy 
is  most  likely  due  to  the  different  types  of  forcing  mechanisms  used  in  these  experiments. 
Both  the  investigations  of  Strykowski  and  Sreenivasan  and  of  Williams  et  al.  employ 
forcing  mechanisms  which  were  very  two-dimensional  in  nature.  The  vibrating  ribbon, 
however,  necessarily  has  a  spanwise  variation  in  forcing  amplitude  of  100  percent  from 
the  maximum  in  the  center  to  zero  at  the  fixed  ends.  This  causes  no  problem,  of  course, 
if  the  resulting  disturbances  are  observed  very  close  to  the  ribbon  and  thus  far  from  the 
region  of  influence  of  the  clamped  ends.  This  was  precisely  the  case  in  the  experiments  of 
Schubauer  and  Scramstad.  If  the  disturbances  are  observed  far  from  the  ribbon,  however, 
the  influence  of  the  finite  ends  cannot  be  ignored. 

A  second  and  equally  important  reason  for  looking  at  the  non-linear  development  of 
two-dimensional  waves  is  that  it  provides  a  point  of  departure  for  the  results  of  the  next 
few  chapters  on  oblique  waves.  As  will  be  seen,  oblique  waves  of  sufficiently  large  ampli- 
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tude  show  interesting  behavior  arising  from  both  the  non-linear  nature  of  the  problem  as 
weli  as  the  three  dimensional  geometry  of  the  waves.  Having  first  studied  the  non-linear 
effects  in  the  two  dimensional  case,  it  will  be  much  easier  to  recognize  and  identify  the 
role  of  these  two  distinct  mechanisms  in  the  oblique  wave  development. 

3.2  Experimental  Conditions 

In  the  present  experiments,  much  care  was  taken  to  insure  that  the  forced  distur¬ 
bances  were  as  two-dimensional  as  possible.  Two  sources  of  three-dimensionality  were 
inherent  in  the  forcing  mechanism  employed  in  these  experiments.  As  pointed  out  previ¬ 
ously,  the  finite  segmentation  of  the  heater  array  introduces  higher  spanwise  wavenumber 
components  into  the  initial  disturbance  waveform.  These  are  expected  to  be  strongly 
damped  in  accordance  with  the  linear  theory,  and  the  departure  from  two-dimensionality 
due  to  this  effect  should  therefore  be  small  far  downstream  of  the  heater. 

A  second  source  of  three-dimensionality  is  the  variation  in  the  resistance  of  the 
heaters  themselves.  Though  great  care  was  taken  in  the  array  fabrication,  the  individ¬ 
ual  heater  resistances  can  vary  by  as  much  as  ten  percent  from  the  mean.  The  heat 
flux  delivered  to  the  flow,  however,  can  be  made  constant  across  the  span  of  the  array 
by  adjusting  the  forcing  voltage  via  a  variable  resistor  (one  per  heater)  in  the  control 
electronics.  In  this  way,  the  power  delivered  to  the  flow  could  be  maintained  at  a  very 
nearly  constant  level  across  the  35  cm  span  of  the  tunnel. 

3.3  Experimental  Results  for  Single  Mode  Excitation 
3.3.1  Disturbance  Growth  Rates 

Figure  3.1  shows  a  typical  wall  shear  time  series  of  the  naturally  occurring  T.  S. 


waves.  In  the  upper  plot,  the  fluctuating  wall-shear  r'  is  plotted  normalized  by  the 
mean  wall  shear  f  and  has  an  rms  value  of  less  than  one  percent.  Plotted  below  is 
the  amplitude  spectrum.  Nearly  all  of  the  energy  is  found  between  20  and  30  hz  (non- 
dimensional  frequency,  F  =  0.78  to  1.18).  Virtually  no  energy  is  seen  outside  of  this 
narrow  frequency  band.  With  the  forcing  turned  on,  the  observed  wall  shear  is  as  shown 
in  figure  3.2.  In  this  particular  case  the  forcing  frequency  is  22  hz  (F  =  0.86),  and  the 
power  supplied  to  the  32  element  heater  array  is  approximately  460  watts  (18  watts 
per  square  cm).  The  forced  waves  are  seen  to  be  quite  regular,  with  the  variation  in 
amplitude  due  almost  entirely  to  the  naturally  occurring  background  waves. 

In  order  to  verify  the  two-dimensionality  of  the  forced  oscillations,  two  hot-film 
wall-shear  sensors  at  spanwise  locations  z  =  ±2 .0cm  were  used.  As  seen  in  figure  3.3, 
the  waves  are  indeed  quite  two-dimensional.  There  is  a  slight  difference  in  disturbance 
amplitude  at  the  two  probes,  but  it  is  generally  less  than  a  few  percent  of  the  mean 
value.  One  possible  reason  for  this  difference  would  be  a  variation  in  the  amplitude  of 
the  naturally  occurring  waves  at  each  sensor  location.  Also,  a  slight  mismatch  in  the 
trimming  resistors  in  the  control  electronics  could  account  for  this  difference. 

In  order  to  assess  the  extent  of  non-linearity  in  the  disturbance  development,  it  is 
instructive  to  look  at  the  wave  amplitude  as  a  function  of  the  power  input  to  the  heater 
array.  Figure  3.4  shows  the  normalized  disturbance  amplitude  at  two  sensor  locations 
as  a  function  of  the  input  power.  It  is  noted  that  the  disturbance  amplitude  is  a  very 
nearly  linear  function  of  the  heat  input  up  to  forcing  levels  of  about  300  W.  Beyond  this 
forcing  level,  a  departure  from  linearity  is  observed.  The  linearity  of  the  fluid  response 
is  in  agreement  with  equation  1.16  provided  that  the  quantity  dp/dT  remains  constant. 
At  the  higher  forcing  levels,  however,  the  tomperature  difference  AT  as  obtained  from 
equation  1.11  becomes  quite  large  (AT  =  20°C  for  Qt  =  500W),  and  the  change  in 
dft/dT  over  this  temperature  range  can  no  longer  be  neglected.  For  a  temperature 
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increase  of  20°C,  dft/dT  decreases  by  nearly  fifty  percent.  This  change  in  the  fluid 
properties  with  temperature  therefore  accounts  for  the  change  in  the  magnitude  of  the 
velocity  perturbation  introduced  into  the  fluid. 

It  is  important  at  this  point  to  differentiate  between  the  mechanism  by  which  the 
disturbances  are  introduced  into  the  fluid  and  the  mechanism  of  viscous  instability  which 
governs  their  subsequent  development.  As  has  just  been  seen,  the  velocity  perturbation 
is  introduced  into  the  fluid  by  a  non-linear  mechanism.  Even  at  the  largest  forcing 
amplitudes,  however,  the  velocity  perturbation  which  is  initially  introduced  into  the 
fluid  is  very  small.  The  downstream  development  of  the  disturbance,  then,  is  governed 
by  the  linear  theory.  It  is  the  effect  of  weak  non-linearity  in  this  region  which  is  of 
primary  interest. 

In  the  absence  of  non-linear  effects,  the  disturbance  amplitude  measured  at  a  given 
location,  z,  would  be  given  by 

A(x)  —  A(*o)f  e~aiXdx  (3.1) 

J  *u 

where  a,  is  the  spatial  growth  rate  given  by  the  linear  theory,  and  A(zo)  is  a  reference 
amplitude  measured  at  z<j.  Thus,  if  we  examine  the  amplitudes  at  two  different  sensors 


located  at  streamwise  positions  zo  and  Zi,  the  amplitude  ratio  A(x\)/A(xo)  should  be 
a  constant  value.  Any  departure  from  this  value  would  indicate  a  departure  from  the 
linear  theory.  Figure  3.5  shows  a  plot  of  such  an  amplitude  ratio  evaluated  at  Reynolds 
numbers  Ref  —  1178  and  1240.  The  amplitude  ratio  A\/Aq  is  plotted  as  a  function  of 
the  amplitude  Ai  evaluated  at  Rtf  =  1240.  As  can  be  seen  in  figure  3.5,  the  amplitude 
ratio  is  quite  constant  up  to  disturbance  amplitudes  of  nearly  12  percent.  Substantial 
variation  is  seen  at  the  lower  disturbance  amplitudes  since  the  forced  oscillations  are 
on  the  order  of  those  naturally  ocurring,  and  the  resulting  growth  rate  depends  on  the 
exact  phase  relation  between  these  two.  As  the  amplitude  of  the  forced  oscillations 
increases,  however,  the  contribution  of  the  background  waves  becomes  less  important. 


It  is  thus  concluded  that  the  growth  rates  as  predicted  by  the  linear  theory  remain  valid 
for  disturbance  amplitudes  of  at  least  12  percent. 


3.3.2  Generation  of  Higher  Harmonica 

Returning  to  the  spectrum  of  figure  3.2,  one  further  feature  is  observed.  A  small  but 
clearly  evident  peak  exists  at  44  hz  (F=1.72),  exactly  twice  the  forcing  frequency.  This 
peak  was  not  apparent  in  the  unforced  case  (figure  3.1).  The  nature  of  this  component  of 
the  disturbance  can  be  seen  by  narrowly  band-pass  filtering  the  time  series  data  around 
the  two  dominant  spectral  peaks  at  22  Hz  and  44  Hz  and  superposing  these  two  filtered 
time  series  traces  as  shown  in  figure  3.6.  The  smaller  amplitude  44  Hz  trace  can  now  be 
identified  as  a  non-linearly  generated  first  harmonic  of  the  fundamental  T.  S.  wave  at  22 
Hz.  This  is  evident  by  noting  that  the  two  time  series  traces  are  perfectly  phase-locked. 
If  the  wave  at  44  Hz  were  a  linear  wave  itself,  the  wavespeed  would  be  different  from 
that  of  the  22  Hz  wave  and  phase-locking  would  not  occur. 

Further  insight  into  the  nature  of  the  first  harmonic  is  obtained  by  examinating  its 
amplitude.  This  is  shown  in  figure  3.7  as  a  function  of  the  forcing  level.  The  amplitude 
is  generally  seen  to  increase  with  the  forcing  level,  but  the  important  point  to  notice  here 
is  that  the  amplitude  is  increasing  with  the  Reynolds  number  as  well.  According  to  the 
linear  theory,  a  disturbance  at  this  frequency  and  Reynolds  number  should  be  strongly 
damped,  i.e.  decreasing  with  increasing  Reynolds  number.  This  is  clearly  shown  in  the 
neutral  curve  of  figure  1.1. 

The  behavior  of  the  first  harmonic  wave  is  therefore  not  described  by  the  linear 
theory.  This  is  not  suprising  since  this  is  in  fact  a  non-linear  phenomenon.  It  arises 
due  to  the  self-interaction  of  the  fundamental  wave  with  itself.  The  amplitude  of  the 
first  harmonic  is  directly  related  to  that  of  the  fundamental  as  seen  in  figure  3.8.  The 


dependence  is  quadratic  as  expected.  Higher  harmonics  are  aiso  generated  but  are  only 
seen  in  this  investigation  when  the  fundamental  amplitude  becomes  very  large.  Figure 

3.9  shows  one  example  of  the  appearance  of  higher  harmonics.  The  fundamental,  first 
harmonic,  and  second  harmonic  are  clearly  visible,  and  a  very  small  peak  at  the  third 
harmonic  is  also  seen. 

3.3.3  Frequency  Dependence 

In  the  previous  section,  it  was  shown  that  the  integrated  amplification  A\/Aq  was 
independent  of  the  disturbance  amplitude  at  least  up  to  amplitudes  of  twelve  percent. 
This  is  a  local  result  in  that  we  have  examined  the  wave  behavior  over  a  very  narrow 
range  of  Reynolds  numbers.  Furthermore,  we  have  so  far  studied  disturbances  at  only 
one  frequency. 

In  order  to  obtain  a  more  global  view,  the  frequency  dependence  of  the  forced 
oscillations  was  studied  by  performing  the  following  experiment.  The  phased  heater 
array  was  initially  programmed  to  force  two-dimensional  waves  of  small  amplitude.  The 
frequency  of  the  forced  disturbance  was  then  increased  in  increments  of  2  Hz,  completely 
covering  the  range  of  unstable  frequencies  seen  by  a  sensor  at  Res-  «  1240.  In  this  way, 
the  integrated  disturbance  amplitude  is  obtained  as  a  function  of  frequency.  The  output 
voltage  of  the  heater  array  was  then  increased,  and  the  experiment  was  repeated.  Figure 

3.10  shows  the  result  of  this  experiment,  where  the  normalized  fluctuating  wall-shear  t' /f 
is  plotted  as  a  function  of  both  the  frequency  and  the  forcing  level.  At  each  forcing  level, 
a  Gaussian  least  squares  fit  is  shown  as  well.  There  is  no  a  priori  reason  for  choosing  a 
Gaussian  other  than  that  it  seems  to  fit  the  data  reasonably  well. 

The  important  point  to  be  noticed  from  figure  3.10  is  that  the  region  of  unstable 
frequencies  as  seen  by  the  wall-shear  sensor  is  virtually  independent  of  the  amplitude  of 
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the  forced  waves.  At  each  forcing  level,  the  most  amplified  frequency  is  25  hz,  and  the 
range  of  unstable  frequencies  extends  from  12  hz  to  40  hz.  Since  figure  3.10  gives  the 
integrated  amplification  from  the  heater  to  the  sensor,  it  seems  reasonable  to  conclude 
that  the  growth  rates,  a,,  are  independent  of  the  disturbance  amplitude  for  a  very  wide 
range  of  frequencies  and  Reynolds  numbers.  Thus  if  one  were  to  inquire  as  to  the  shape 
of  the  neutral  curve  for  finite  amplitude  waves,  the  present  results  suggest  that  it  is 
essentially  unchanged  up  to  disturbance  amplitudes  of  12  percent. 

3.3.4  Summary  of  Observed  Non-linear  Effects 

The  effect  of  weak  non-linearity  in  boundary  layer  transition  is  thus  seen  to  be 
relatively  benign.  The  growth  rates  were  seen  to  follow  the  linear  theory  up  to  pertur¬ 
bation  amplitudes  of  nearly  twelve  percent.  This  was  seen  by  analyzing  both  the  local 
growth  rates  as  well  as  the  integrated  amplification  over  a  wide  range  of  frequencies  and 
Reynolds  numbers.  The  only  deviation  from  the  linear  theory  arose  in  the  form  of  a 
non- linearly  generated  first  harmonic  wave  phase-locked  to  the  fundamental,  whose  am¬ 
plitude  depended  quadratically  upon  that  of  the  fundamental.  Higher  harmonics  were 
also  observed,  but  as  their  amplitudes  scale  with  higher  powers  of  the  fundamental  ampli¬ 
tude,  they  are  only  seen  at  very  large  forcing  levels.  The  appearance  of  these  non-linear 
harmonics  will  certainly  alter  the  wave  behavior  somewhat,  but  as  will  be  seen  in  the 
next  few  chapters,  the  onset  of  three-dimensionality  has  a  much  more  profound  impact 
on  the  transition  process. 


CHAPTER  4 


GEOMETRICAL  NATURE  OF  OBLIQUE  INSTABILITY  WAVES 
4.1  Introduction 

Before  we  begin  to  examine  the  behavior  of  oblique  instability  waves  in  a  boundary 
layer  flow,  it  is  necessary  that  we  first  define  the  type  of  disturbance  that  is  to  be  studied. 
Throughout  the  next  several  chapters,  the  disturbance  velocity  field  will  be  assumed  to 
be  of  the  form 

G(x,t)  =  u(y)e,(a*+**-“'‘)  (4.1) 

The  velocity  components  u(x,  t)  are  complex;  all  experimental  data  will  of  course  corre¬ 
spond  to  only  the  real  part,  however.  The  angular  frequency  w  =  2nf  will  be  taken  as 
a  purely  real  quantity  corresponding  to  the  physically  realistic  spatial  growth  problem. 
Most  analytical  studies,  on  the  other  hand,  consider  the  temporal  growth  problem  where 
oj  is  complex.  The  streamwise  wavenumber  a  is  complex  with  a,  giving  the  rate  of  growth 
(a,  <  0)  or  decay  (a<  >  0)  of  the  disturbance  amplitude.  The  real  part  ar  =  2*7 A* 
is  the  streamwise  wavenumber  of  the  disturbance.  The  spanwise  wavenumber  0  is  in 
general  also  complex,  but  under  the  present  experimental  conditions  it  can  be  shown  to 
be  purely  real;  i.e.  0  =  0r  =  2x/A,.  This  will  be  shown  in  the  next  section. 

Equation  4. 1  is  the  most  commonly  used  form  for  oblique  waves,  but  other  notations 
are  frequently  also  found  in  the  literature.  For  example,  equation  4.1  could  equally  well 
be  written  as 

u(x,t)  =  (4.2) 

a  a  a  a 

where  Z  =  at  +  0k  a  the  total  wavenumber  vector,  and  c  =  cxi  +  c,k  is  the  phase 
velocity  vector  of  the  disturbance.  Written  in  this  form,  one  must  ask  the  question:  In 
what  direction  does  the  disturbance  travel?  That  is,  we  seek  the  components  of  the  phase 
velocity  vector,  c*  and  cM.  For  that  matter,  one  would  also  like  to  know  the  magnitude 
and  direction  of  the  group  velocity  vector  C}  =  dw/dn,  as  well,  as  this  quantity  gives 


the  direction  in  which  energy  is  propagated.  These  quantities  are  of  little  importance  in 
theoretical  studies,  since  the  disturbances  are  usually  assumed  infinitely  periodic  in  the 
x  and  z  directions  and  are  thus  prevented  from  “propagating”  out  of  the  computational 
domain.  In  an  experiment,  however,  the  direction  of  propagation  must  be  known  in 
order  to  correctly  position  sensors  for  measurement.  Determination  of  the  group  and 
phase  velocities  is  thus  an  important  practical  matter,  and  it  will  be  shown  in  section 
4.3  that  an  understanding  of  these  quantities  adds  a  good  deal  of  physical  insight  into 
the  nature  of  oblique  instability  waves  as  well. 

Finally,  section  4.4  will  take  a  look  at  the  vorticity  field  of  a  single  oblique  mode. 
It  will  be  shown  that  whereas  the  vorticity  field  of  a  two-dimensional  disturbance  is 
purely  one-dimensional  (i.e.  a  scalar,  ojm),  the  vorticity  field  for  an  oblique  normal  mode 
is  fully  three-dimensional  with  <*/a,  <«/y,  and  ujb  each  being  a  non-vanishing  quantity. 
The  implications  of  this  revealed  three-dimensional  structure  on  the  stability  of  the 
disturbances  will  be  discussed  as  well. 


4.2  Spatial  Growth  Rates 


In  the  words  of  Mack  (1984),  "It  is  also  possible  to  conceive  of  wavemakers  that 
excite  single  oblique  normal  modes  in  boundary  layers  that  are  independent  of  z.  Such 
normal  modes  will  have  an  initial  0T  which  matches  that  of  the  wavemaker,  and,  because 
the  wave  can  grow  only  in  x,  the  initial  0%  must  be  zero.”  The  phased  heater  array  used 
in  the  present  study  is  precisely  such  a  wavemaker. 


It  is  easy  to  see  why  0*  must  be  identically  zero  by  performing  the  following  thought 
experiment.  Imagine  a  heater  array  such  as  the  one  described  in  Chapter  2  but  with 
a  span  extending  to  infinity  in  both  directions.  Assume  now  that  the  elements  of  such 
an  array  are  excited  with  the  same  amplitude  and  phased  so  as  to  produce  a  single 


oblique  normal  mode.  A  sensor  located  some  distance  downstream  of  this  heater  would 
measure  the  same  r.m.s.  disturbance  amplitude  everywhere  across  the  span  as  there  is  no 
preferred  spanwise  direction  and  no  definable  spanwise  origin  for  the  disturbance.  The 
present  heater  array  completely  spans  the  test  facility  very  nearly  approximating  such 
a  hypothetical  array.  Thus  for  the  present  experimental  conditions,  we  may  justifiably 
assume  that  /3,  =  0. 


4.3  On  the  Direction  of  Propagation 


The  group  and  phase  velocities  are  important  quantities  in  any  wave-like  phe¬ 
nomenon.  An  examination  of  the  current  literature  shows  that  there  is  some  confusion 
regarding  these  quantities.  For  example,  Craik  (1971)  considers  a  disturbance  of  the 
form  4.1,  with  w  complex  and  a  real,  and  assumes  that 


That  is,  the  phase  velocity  is  assumed  to  be  parallel  to  the  freestream  direction.  Mack 
(1984),  on  the  other  hand,  considers  the  very  same  disturbance  and  assumes  that 


(4.4) 


i.e.  the  phase  velocity  is  assumed  to  be  parallel  to  the  wavenumber  vector  k.  Similarly, 
Kachanov  and  Levchenko  (1983)  take  c  parallel  to  k  in  the  analysis  of  their  experimental 
results. 


Ip  an  attempt  to  resolve  this  issue,  the  following  experiment  was  conducted.  A  single 
oblique  wave  of  finite  span  was  generated  by  the  phased  heater  array.  The  span  of  the 
disturbance  was  35  cm,  i.e.  the  full  span  of  the  test  section.  To  study  the  effects  of  finite 
span  as  well  as  the  direction  of  propagation,  the  disturbance  was  then  “shifted”  across 
the  span  of  the  array.  The  finite  disturbance  was  “shifted”  by  turning  heater  elements 


-  33- 

off  in  a  sequential  manner  beginning  at  one  end  of  the  array  and  proceeding  across  the 
array  to  the  other  end.  As  far  as  the  sensor  on  the  plate  centerline  is  concerned,  this 
procedure  is  the  same  as  shifting  the  disturbance  through  the  tunnel  side  wall,  since 
disturbances  near  the  side  walls  have  no  effect  on  the  flow  at  the  plate  centerline.  In  this 
way  a  single  probe  located  on  the  plate  centerline  was  used  to  record  the  entire  spanwise 
extent  of  the  disturbance.  This  procedure  was  then  repeated,  “shifting”  the  disturbance 
in  the  opposite  direction;  i.e.  “through”  the  other  tunnel  side  wall. 

The  spanwise  amplitude  distribution  obtained  in  this  manner  gives  results  which 
are  comparable  to  those  that  would  be  obtained  by  moving  a  sensor  relative  to  a  fixed 
disturbance  of  essentially  twice  the  span.  Thus,  we  have  effectively  doubled  the  width 
of  the  tunnel.  The  advantages  of  this  scheme  are  twofold.  First,  the  edge  contamination 
region  near  the  tunnel  side  walls  is  completely  avoided.  Similarly,  the  finite  end  effects  of 
the  disturbance  itself  are  kept  to  a  minimum  by  making  the  disturbance  width  as  large 
as  possible.  This  leaves  a  substantial  segment  in  the  center  of  the  disturbance  which 
is  not  affected  by  the  finite  ends.  It  is  the  disturbance  in  this  region  that  is  of  most 
interest,  as  it  corresponds  very  nearly  to  a  single  oblique  normal  mode. 


Figure  4.1a  shows  the  amplitude  distribution  obtained  in  this  manner.  Several  ob¬ 
servations  can  be  made  from  this  plot.  The  disturbance  is  seen  to  extend  approximately 
20  cm  to  either  side  of  the  plate  centerline  with  a  relatively  constant  amplitude  section 
from  -10  cm  to  10  cm.  A  significant  amount  of  scatter  exists  in  the  amplitude  and  is  ac¬ 
counted  for  by  the  naturally  occurring  background  oscillations.  If  we  assume  a  spreading 
angle  of  16  degrees  for  the  finite  end  effects  as  suggested  by  Mack  (1984),  then  spanwise 
locations  where  |s|  >  20  should  be  free  of  the  influence  of  the  forced  oscillations.  The 
scatter  in  this  region  of  the  flow  thus  represents  the  scatter  in  the  naturally  occurring 
oscillations. 


The  angle  of  the  oblique  wave  in  this  experiment  was  approximately  15  degree: 


I, 


'« 


The  streamwise  distance  that  the  disturbance  has  travelled  from  its  origin  at  the  heater 
to  the  point  of  measurement  is  over  30  cm.  From  figure  4.1a,  it  is  seen  that  the  energy 
of  the  disturbance  has  travelled  very  nearly  in  the  stream  direction,  spreading  perhaps 
a  few  centimeters  further  to  values  of  positive  z  than  to  values  of  negative  z.  This 
slight  asymmetry  will  be  discussed  later,  but  for  the  present,  it  seems  reasonable  to 
conclude  that  the  group  velocity  of  this  disturbance,  or  the  velocity  with  which  energy  is 
propagated,  is  very  nearly  oriented  in  the  freestream  direction.  It  is  certainly  not  parallel 
to  the  wavenumber  vector  k,  as  propagation  in  that  direction  would  have  displaced  the 
energy  distribution  nearly  9  cm  from  the  plate  centerline.  Furthermore,  since  the  center 
region  of  the  disturbance  consists  of  a  single  normal  mode,  the  phase  velocity  in  this 
region  must  be  equal  to  the  group  velocity.  This  must  be  the  case,  since  there  can  be  no 
possible  dispersion  for  a  single  mode.  We  therefore  conclude  that  the  phase  velocity  is 
oriented  in  the  freestream  direction  as  well. 

An  interesting  comparison  with  the  preceding  experiment  is  provided  by  the  point 
source  calculations  of  Mack  (1984).  Figure  4.1b  shows  a  plot  of  the  normalized  dis¬ 
turbance  amplitude  of  a  single  oblique  normal  mode  of  finite  span.  The  origin  of  the 
disturbance  in  this  case  was  physically  oblique;  i.e.  one  end  of  the  disturbance  was  lo¬ 
cated  at  a  higher  Reynolds  number  than  the  other  end.  As  in  figure  4.1a,  a  region  in 
the  center  exists  which  is  free  of  the  effects  of  the  finite  ends.  In  this  case,  however,  the 


amplitude  is  not  constant  but  rather  increases  monotonically  across  the  span.  This  is  due 
to  the  physically  oblique  nature  of  the  disturbance  origin.  The  disturbance  has  simply 
travelled  diferrent  distances  and  therefore  undergone  a  different  integrated  amplification 
at  each  spanwise  location.  It  is  also  to  be  noted  that  the  disturbance  amplitude  falls  off 
rather  more  smoothly  on  the  left  side  of  each  figure  and  is  seen  to  fluctuate  somewhat 
more  on  the  right  side.  Surprisingly,  this  is  a  very  repeatable  feature  of  the  experimental 
data. 
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A  plausible  explanation  for  this  feature  is  offered  by  noting  that  a  sum  over  oblique 
normal  modes  as  in  equation  1.18  would  have  to  be  constructed  in  these  regions  in 
order  to  properly  represent  the  finite  nature  of  the  disturbance  span.  As  will  be  seen 
later  in  this  chapter,  the  vortex  filaments  of  such  a  disturbance  are  oriented  very  nearly 
perpendicular  to  the  wavenumber  vector  and  thus  lay  at  an  oblique  angle  to  the  flow 
direction.  By  analogy  with  the  starting  vortex  of  a  finite  span  airfoil,  these  vortices 
must  also  connect  to  some  form  of  starting  vortex  system  at  the  heater.  The  more 
downstream  end  of  the  oblique  vortices  must  therefore  turn  through  a  greater  angle 
in  order  to  connect  with  a  starting  vortex  system  than  the  more  upstream  end.  Thus 
the  sum  of  equation  1.18  would  consist  of  different  combinations  of  spanwise  modes  for 
the  two  finite  ends  of  the  disturbance.  The  more  downstream  end  would  require  more 
energy  at  higher  spanwise  wavenumbers  in  order  to  turn  the  corner  and  could  thereby 
account  for  the  observed  short  wavelength  oscillation  on  the  right  side  of  both  figure 
4.1(a)  and  4.1(b)  as  well  as  the  slight  difference  in  lateral  spreading  of  the  disturbance. 
It  must  be  remembered,  however,  that  this  argument  is  strictly  valid  only  for  an  inviscid 
fluid,  whereas  the  present  flow  is  viscous.  It  is  nevertheless  believed  to  provide  a  useful 
explanation  for  the  finite  span  vorticity  field  in  the  present  experiments. 


4.3.1  The  Nature  of  Viscous  Instability 

In  the  previous  section,  the  direction  of  propagation  of  an  oblique  wave  was  studied. 
It  is  also  of  interest  to  examine  the  magnitude  of  this  propagation  velocity.  In  particular, 
it  is  interesting  to  see  how  this  quantity  varies  with  the  angle  of  the  wave.  Consider  the 
following  experiment.  The  phased  heater  array  was  initially  programmed  to  force  two- 
dimensional  waves.  Two  hot-film  wall-shear  sensors  separated  by  3  mm  in  the  streamwise 
direction  were  used  to  obtain  the  phase  velocity.  The  time  of  flight  of  individual  wave 


crests  was  measured  and  averaged  over  50  cycles  of  the  disturbance.  Dividing  the  sensor 
separation  distance  by  the  time  of  flight  gives  the  phase  speed  of  the  disturbance.  The 
heater  array  was  then  programmed  to  output  oblique  waves  at  increasingly  greater  angles 
to  the  flow  direction,  and  each  time  the  phase  speed  was  measured. 

The  result  of  this  experiment  is  shown  in  figure  4.2.  The  normalized  phase  velocity 
Cx/Uoa  is  plotted  as  a  function  of  the  spanwise  wavenumber  /?.  Surprisingly,  the  phase 
velocity  is  seen  to  remain  constant.  That  is,  the  propagation  velocity  cx  is  independent 
of  the  angle  of  the  disturbance.  In  order  to  see  why  this  should  be  the  case,  let  us 
examine  the  physical  mechanism  governing  viscous  instability. 

The  physical  mechanism  governing  these  instability  waves  can  be  viewed  as  follows. 
Sinusoidal  oscillations  in  the  velocity  profile  cause  vorticity  to  be  generated  at  the  wall. 
This  vorticity  then  diffuses  away  from  the  wall  and  is  convected  by  the  mean  velocity 
field  in  the  streamwise  direction.  Convection  by  the  disturbance  velocity  field  can  be 
neglected,  as  it  is  zero  in  the  m^an.  When  the  vorticity  from  the  wall  has  diffused  to 
the  critical  layer,  it  adds  to  the  vorticity  which  is  already  present  in  the  disturbance.  If 
the  frequency  of  oscillation  and  critical  layer  height  yer(Re)  are  such  that  the  vorticity 
shed  from  the  wall  arrives  at  yer  in  phase  with  that  already  present,  then  the  vorticity  of 
the  disturbance  increases.  If,  on  the  other  hand,  the  shed  vorticity  is  out  of  phase  upon 
reaching  the  critical  layer,  then  the  vorticity  of  the  disturbance  decreases.  Viewed  in  this 
way,  it  is  seen  that  viscous  instability  waves  are  in  a  sense  not  waves  at  all.  It  was  seen 
earlier  that  these  waves  do  not  propagate  in  the  direction  normal  to  their  phase  fronts 
as  normal  waves  do,  but  rather  travel  in  the  freestream  flow  direction.  It  is  perhaps  not 
even  correct  to  speak  of  wave  “propagation”  at  all,  since  the  problem  is  really  one  of 
diffusion  and  convection  of  periodic  vorticity  rather  than  propagation  of  a  wave. 

When  the  vortices  are  oblique,  the  situation  is  basically  the  same.  Now,  oscillations 
of  the  velocity  profile  cause  both  x-vorticity  and  z- vorticity  to  be  created  at  the  wall 
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This  vorticity  again  diffuses  away  from  the  wall  and  is  convected  by  the  mean  flow  in 
the  x-direction.  Therefore,  it  is  not  really  suprising  that  the  phase  speed  of  the  waves 
is  seen  to  be  independent  of  the  wave  angle  as  in  figure  4.2,  for  the  same  mechanisms  of 
diffusion  and  convection  of  vorticity  govern  the  disturbance  development  regardless  of 
the  orientation  of  the  lines  of  constant  phase. ^ 


4.4  A  Look  at  the  Vorticity  Field 

A  great  deal  of  insight  can  be  gained  by  looking  at  the  vorticity  field  of  the  pertur¬ 
bation.  For  example,  in  the  two  dimensional  case,  the  vorticity  field  is  seen  to  consist  of 
a  streamwise  periodic  array  of  vortices  with  the  only  non-vanishing  vorticity  component 
being  that  in  the  spanwise  direction, 

dV  dU  I'  dU\  i(ai-ut)  (A 

"* =  ai  -  ai  =  <K"  ‘  (4'5> 

In  the  case  of  oblique  disturbances  of  the  form  4.1,  it  is  natural  to  expect  that  the 
vorticity  field  would  be  very  similar  to  that  of  a  two-dimensional  disturbance,  with  the 
only  difference  being  that  the  vortex  filaments  are  rotated  by  an  angle  0  =  tan_,(/?/a). 
Thus,  we  would  expect  two  vorticity  components,  u>M  and  wx.  This  picture  of  the  vorticity 
field  of  an  oblique  wave  is,  however,  not  entirely  correct.  To  see  why  this  physical  picture 
is  incorrect,  it  is  necessary  to  look  at  the  mathematical  form  of  the  disturbance  velocity 
field  as  given  in  equation  4.1.  Taking  the  curl  of  4.1,  the  vorticity  field  is  obtained  as 


“■  -  %  -  T.  ■  -  tf******— ' 

w,  =  ^  -  ~~  =  (i/?ti  -  ia 

«.-£-£  =  (ft*  - 

ox  dy  '  dy' 


(4.6a) 


(4.66) 


(4.6c) 


S 


The  difference  between  the  two-dimensional  and  three-dimensional  vorticity  fields 
now  becomes  apparent.  In  our  previous  picture  of  the  three-dimensional  vorticity  field, 


the  vortex  filaments  were  straight  and  everywhere  perpendicular  to  the  wavenumber 


vector  k  =  ai+0k.  More  importantly  though,  this  picture  of  the  vorticity  field  contained 


no  vertical  vorticity  uv,  whereas  from  4.6b,  it  is  seen  that  wy  is  not  identically  zero. 


To  prove  that  u>y  is  in  fact  a  non-vanishing  quantity  for  an  oblique  wave,  consider 


the  linearized  equations  of  motion: 


V  •  u  =  0 


(4.7a) 


,d  „  d  w  dl/,  _  1  _o  — 

^+^)U  +  W  =  _VP+^VU 


(4.76) 


Substituting  the  assumed  form  of  the  disturbance  from  4.1  and  correspondingly  taking 


the  pressure  perturbation  as 


p(*,0  =  p(yK(~+*'-u’t) 


we  obtain 


dv 

tau  +  —  +  ipw  =  0 
dy 


(4.9a) 


r  •  .  „  1  .  d2  -  dU 


(4.96) 
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(4.9c) 


[-tw  +  iaU  -  -  a2  -  0*)]w  =  -*Pp 


(4.9d) 


Multiplying  4.9(b)  by  */0,  4.9(d)  by  ia,  and  subtracting,  we  obtain 


[»w  -  iaU  +  -  <*2  "  £2)](«7?«  ~  *aw)  = 


(4.10) 


The  fact  that  uv  is  not  identically  zero  then  follows  by  contradiction,  for  assuming  that 


wy  =  0  implies  that  the  (eft  hand  side  of  4.10  vanishes.  Since  dU/dy  £  0  and  #  0, 


this  implies  that  v  must  be  zero  everywhere.  This  problem  was  considered  in  1933  by 


Squire,  where  he  showed  that  such  a  disturbance  (an  oblique  wave  with  v(y)  =  0)  is 
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always  stable.  Therefore  if  uv  =  0,  no  instability  exists  for  oblique  waves  at  all.  Since  it 


is  known  from  experiments  that  this  is  not  the  case,  it  must  be  concluded  that  u >y  is  not 


identically  zero  and  is  in  fact  given  by  the  expression  4.6(b). 


4.5  Three-Dimensional  Stability  Theory  Revisited 


Having  determined  that  the  vorticity  field  of  an  oblique  wave  is  indeed  fully  three- 


dimensional,  let  us  now  examine  the  possible  impact  that  this  might  have  on  the  stability 


of  the  disturbance.  It  is  convenient  to  begin  with  the  full  (non-linearized)  equations  of 


motion  expressed  in  terms  of  the  vorticity: 


—  +  (u  •  V)<3  =  (w  •  V)u  +  i/V2w 
at 


(4.11) 


Here  (2  represents  the  total  vorticity  of  the  flow,  i.e.  the  mean  vorticity  plus  that  of 


the  disturbance.  The  velocity  vector  U  is  likewise  that  of  the  mean  flow  plus  that  of 


the  disturbance.  Each  term  in  4.11  is  easily  interpreted  as  follows:  dJ/dt  is  the  rate  of 


change  of  vorticity,  (tx  •  V)w  is  the  convection  of  vorticity  by  the  velocity  field,  (w  •  V)« 


corresponds  to  the  rotation  and  stretching  of  vorticity  by  the  strain  field,  and  the  final 


term  is  of  course  the  dissipation  due  to  viscosity.  It  is  the  first  term  on  the  right-hand 


side  that  is  of  particular  interest  here,  as  it  was  identically  zero  in  the  two-dimensional 


case.  For  oblique  waves,  however,  it  provides  an  alternate  means  of  vorticity  production 


that  was  not  available  for  two-dimensional  disturbances. 


It  is  easily  shown  that  equation  4.11  can  be  rewritten  as 


—  +  (u  ■  V)w  =  <J  •  e  +  i/V2uj 
at 


(4.12) 


where  the  strain  rate  tensor  t  is  given  by 


1 ,  du,  duj  x 
e,J  2  dx}  +  dxi 


(4.13) 
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If  we  now  linearize  4.12,  making  the  usual  parallel  flow  assumption,  U  =  U(y),  we 


obtain 


d<2'  rtd<2'  ,d*U 
-rr  +  U  — - v 


=  q'  •  i  +  n  •  e'  +  i/V 


V 


(4.14) 


dt  dx  dy3 

where  e  =  ?+«'  and  (2  —  Q+cJ'.  The  vorticity  production  term,  (<J- V)u  =  c5*t,  has  thus 
been  separated  into  the  two  linear  components  shown  on  the  right-hand  side  of  4.14. 
The  first  of  these  corresponds  to  the  distortion  of  the  disturbance  vorticity  field  by  the 
mean  strain  field,  and  the  second  term  gives  the  distortion  of  the  mean  vorticity  field  by 
the  disturbance  strain  field. 


In  order  to  explore  the  nature  of  these  vorticity  production  terms,  let  us  look  at 
the  simplest  possible  problem,  that  of  Couette  flow.  Couette  flow  is  known  to  be  stable 
to  small  two-dimensional  disturbances,  but  the  three-dimensional  problem  remains  rel¬ 
atively  unexplored.  The  mean  flow  profile  is  linear,  and  the  third  term  on  the  left-hand 
side  of  4.14  therefore  vanishes.  If  we  confine  our  attention  for  the  moment  to  only  the 
first  of  the  vorticity  production  terms  and  take  nV  to  be  zero,  we  arrive  at  the  equation 
(dropping  primes  on  the  disturbance  quantities): 

— -  —  w  ■  i  +  i/V2uj  (4.15) 


where  D/ Dt  =  (d/dt  +  Ud/dx).  The  mean  strain  rate  tensor  i  is  given  by 


(4.16) 


The  eigenvalues  of  i  are  found  to  be  dU/dy ,  -dU /dy,  and  0.  The  corresponding 
principal  directions  are  ej  =  (e,  4-  ef ),  =  [ex  -  and  ej  =  et.  Rewriting  equation 
4. 15  in  terms  of  the  vorticity  components  in  the  principal  coordinate  system,  u/t ,  wj ,  and 
u>3,  we  obtain  the  following  set  of  uncoupled  equations  for  the  vorticity  components: 


(4.17c) 


From  equations  4.17,  the  nature  of  the  additional  vorticity  production  terms  now 
becomes  more  apparent.  In  equation  4.17(a),  for  example,  it  can  be  seen  that  following 
the  fluid,  there  are  two  competing  mechanisms  operating.  The  first  term  on  the  right 
hand  side  corresponds  to  the  stretching  of  the  disturbance  vorticity  by  the  mean  strain 
field,  dU/dy,  and  the  second  term  is  the  usual  viscous  dissipation  term.  For  sufficiently 
low  Reynolds  numbers  or  for  dU/dy  — »  0,  the  viscous  term  will  dominate  and  the 
magnitude  of  will  decay  with  time.  On  the  other  hand,  for  large  Reynolds  numbers 
or  large  shear  rates,  the  stretching  term  will  dominate,  and  the  magnitude  of  u>i  will 
increase  exponentially.  The  exponential  growth  would,  of  course,  not  continue  without 
bound,  but  rather  would  give  rise  to  turbulent  flow.  The  qualitative  behavior  of  equations 
4.17  is  therefore  in  agreement  with  the  experimentally  observed  behavior,  for  it  is  known 
that  turbulent  Couette  flow  does  indeed  exist.  In  addition,  in  the  purely  two-dimensional 
case,  only  equation  4.17(c)  remains,  and  all  disturbances  are  damped  following  the  fluid. 
This  is  also  in  agreement  with  the  known  result  that  Couette  flow  is  stable  to  two- 
dimensional  disturbances. 

The  exact  conditions  for  instability  would  of  course  require  that  the  full  linearized 
equation  4.14  be  solved  with  all  terms  included.  The  present  brief  examination  of  the 
disturbance  vorticity  suggests,  however,  that  instability  may  indeed  occur  for  three- 
dimensional  disturbances  even  when  two-dimensional  disturbances  are  stable.  This  is  due 
to  the  additional  vorticity  production  terms  inherent  in  the  three-dimensional  problem. 
In  the  next  chapter,  experimental  evidence  will  be  presented  that  indicates  that  under 


d.a.. - . 
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identical  forcing  conditions,  oblique  waves  can  in  fact  be  much  more  unstable  than  any 


two-dimensional  disturbance. 


Chapter  5 

THREE-DIMENSIONAL,  OBLIQUE  WAVE  EXCITATION 


5.1  Introduction  and  Preliminary  Observations 

In  this  chapter,  we  will  examine  the  results  of  several  experiments  designed  to  ex* 
plore  the  behavior  of  single  oblique  instability  waves.  Such  disturbances  have  been 
discussed  in  the  literature  for  over  50  years  and  are  currently  employed  as  an  integral 
part  of  several  proposed  resonance  mechanisms  as  seen  in  Chapter  1.  Virtually  no  exper¬ 
iments  exist,  however,  which  either  confirm  or  deny  the  predictions  of  the  linear  theory 
with  regard  to  the  behavior  of  oblique  waves.  It  is  the  goal  of  this  section  to  shed  some 
light  on  this  subject.  In  particular,  attention  will  be  focused  on  the  following  issue.  It  is 
the  general  consensus  in  the  literature  that  oblique  waves  are  less  dangerous,  that  is,  less 
unstable  than  their  two-dimensional  counterparts.  The  simple  analysis  of  the  preceding 
chapter,  however,  suggests  that  the  additional  vorticity  production  terms  inherent  in  the 
three-dimensional  problem  may  indeed  result  in  much  higher  growth  rates  than  those 
given  by  the  linear  theory.  If  this  were  to  be  the  case,  then  oblique  waves  should  be  seen 
to  be  more  unstable  than  two-dimensional  waves.  As  will  be  shown,  certain  regimes  of 
the  parameter  space  exist  for  which  this  is  indeed  the  case. 

In  order  to  get  a  feel  for  the  behavior  of  oblique  waves,  consider  the  following  exper¬ 
iment.  The  phased  heater  array  was  programmed  to  force  oblique  waves  at  increasing 
angles  to  the  mean  flow  direction.  The  amplitude  of  the  forcing  was  held  constant 
throughout  as  were  the  freestream  flow  conditions.  At  each  angle,  the  frequency  was 
varied,  completely  spanning  the  range  of  unstable  frequencies.  Figure  5.1  shows  a  plot  of 
the  normalized  disturbance  amplitude  as  recorded  at  Ref  =  1240.  In  order  to  obtain  the 
points  in  figure  5.1,  the  spectra  were  locally  band-pass  filtered  at  the  forcing  frequency. 
The  response  shown  is  thus  the  integrated  flow  response  at  the  forcing  frequency  alone. 


A  flow  response  at  frequencies  other  than  the  forcing  frequency  (the  subharmonic,  for 
example)  would  not  appear  in  figure  5.1.  As  will  be  seen,  the  flow  response  does  indeed 
contain  energy  at  frequencies  other  than  the  forcing  frequency.  This  behavior  constitutes 
a  major  departure  from  the  anticipated  results  and  is  discussed  in  great  detail  in  the  next 
few  sections. 

Several  observations  can  be  made  from  figure  5.1.  First,  it  is  noted  that  the  inte¬ 
grated  flow  response  at  the  forcing  frequency  decreases  with  increasing  wave  angle.  Thus 
it  appears  that  oblique  waves  become  less  unstable  as  the  wave  angle  is  increased.  That 
this  is  not  a  correct  statment  will  be  shown  in  detail  in  the  section  5.3.2.  The  reason  is, 
as  mentioned  above,  that  we  have  not  taken  into  account  the  possible  flow  response  at 
frequencies  other  than  the  forcing  frequency.  For  the  present,  then,  it  can  only  be  stated 
that  the  flow  response  at  the  forcing  frequency  decreases  with  increasing  wave  angle. 

A  second  observation  from  figure  5.1  is  that  approximately  the  same  frequency 
range  is  seen  to  be  unstable  regardless  of  the  wave  angle.  This  suggests  that  the  region 
of  unstable  frequency-Reynolds  number  space  (i.e.  the  neutral  surface)  is  essentially  the 
same  for  all  wave  angles.  It  should  be  recalled  that  in  Chapter  4,  it  was  seen  that  the 
phase  speed  of  oblique  disturbances  was  independent  of  the  wave  angle  as  well.  These 
two  observations  are  consistent  with  the  physical  description  of  viscous  instability  as  a 
balance  between  convection  and  diffusion  of  vorticity.  Thus  it  is  seen  that  the  behavior  of 
oblique  waves  is  quite  similar  to  that  of  two-dimensional  waves.  The  only  difference  seen 
thus  far  has  been  that  the  integrated  amplification  at  the  forcing  frequency  decreases 
with  increasing  wave  angle.  This  is  not  the  only  difference,  however.  As  will  be  seen, 
certain  regions  of  the  parameter  space  exist  in  which  a  different  mechanism  dominates 
the  oblique  wave  development.  This  results  in  a  radical  departure  from  the  type  of 
behavior  seen  in  two-dimensional  waves,  and  is  therefore  the  subject  of  the  remainder  of 
the  chapter. 


5.2  The  Onset  of  Non-stationary  Behavior 


The  remainder  of  this  chapter  breaks  down  into  two  separate  parts.  The  first  half 
is  concerned  with  a  simple  observation  from  the  experimental  data.  That  observation 
is  the  following:  certain  regions  of  the  parameter  space  exist  for  which  oblique  waves 
exhibit  inherently  non-stationary  behavior.  By  non-stationary  we  mean  that  the  spectral 
content  of  the  signal  changes  with  time.  This  behavior  occurs  under  well  controlled 
forcing  conditions  to  be  described  and  is  an  entirely  repeatable  feature  of  the  disturbance. 
Section  5.2.1  will  explore  the  time  series  behavior  in  an  attempt  to  characterize  the 
nature  of  this  non-stationary  behavior.  Section  5.2.2  will  then  discuss  the  subject  of 
spectral  analysis  for  non-stationary  signals,  as  the  usual  methods  such  as  the  Fast  Fourier 
Transform  are  no  longer  applicable.  The  Wigner  distribution  will  be  introduced  as 
a  means  of  non-stationary  frequency  analysis  and  will  be  used  to  support  the  direct 
observations  made  from  the  time  series  data.  Finally,  the  implication  of  this  simple 
observation  on  several  of  the  currently  proposed  resonance  mechanisms  will  be  discussed. 

The  second  half  of  the  chapter  will  present  a  simple  mechanism  which  is  again 
suggested  by  direct  observation  of  the  time  series  data.  This  mechanism  is  simply  that 
the  non-stationary  behavior  is  due  to  the  intermittent  roll-up  or  pairing  of  adjacent 
vortices  of  the  disturbance.  This  mechanism  has  long  been  known  to  play  a  major  role 
in  the  development  of  unbounded  shear  flows  such  as  the  mixing  layer  and  wake.  The 
results  from  a  thorough  experimental  examination  of  the  parameter  space  will  then  be 
discussed  and  will  be  shown  to  be  consistent  with  this  simple  model.  Finally,  the  entire 
picture  of  the  role  of  oblique  waves  in  boundary  layer  transition  will  be  summarized  by 
relating  the  geometrical  properties  revealed  in  chapter  4  to  the  experimental  observations 
and  proposed  mechanism  of  the  present  chapter. 


5.2.1  Time  Series  Behavior 


Figure  5.2  shows  a  typical  wall-shear  time  signal.  The  heater  array  in  this  case  was 
programmed  to  force  a  single  oblique  wave  at  approximately  15  degrees.  The  forcing 
frequency  was  II  Hz,  producing  a  22  Hz  disturbance  due  to  the  squaring  of  the  forcing 
voltage.  The  signal  shown  in  figure  5.2  was  recorded  by  a  hot-film  wall-shear  probe  at 
a  Reynolds  number  of  approximately  1300.  Several  observations  can  be  made  from  this 
plot. 

The  first  point  to  be  noticed  is  that  the  signal  is  quite  non-stationary.  This  is  easily 
seen  by  noting  that  in  the  interval  0.3  <  t  <  0.6  a  single  dominant  frequency  exists  in 
the  wall-shear  fluctuations.  This  is  the  frequency  corresponding  to  the  forcing,  i.e.  22 
hz.  The  r.m.s.  fluctuation  is  typically  on  the  order  of  five  percent  of  the  mean  wall- 
shear.  Moving  to  the  interval  0.6  <  t  <  0.8  a  curious  behavior  is  observed.  Adjacent 
peaks  in  the  time  trace  corresponding  to  local  maxima  in  the  wall-shear  fluctuations 
are  seen  to  merge  yielding  a  signal  which  now  contains  considerable  energy  at  twice  the 
fundamental  period  of  oscillation.  The  result  of  this  pairing  is  further  demonstrated 
in  the  initial  portion  of  the  time  trace,  0.0  <  t  <  0.3,  where  very  large  magnitude 
fluctuations  are  observed  with  a  primary  period  twice  that  of  the  fundamental  period 
of  oscillation.  The  fluctuations  are  seen  to  approach  50  percent  of  the  mean  wall-shear 
level  and  are  much  larger  in  the  positive  direction  than  in  the  negative.  This  asymmetry 
with  respect  to  the  sign  of  the  fluctuations  is  a  very  consistent  feature  of  the  observed 
experimental  data.  These  observations  will  be  shown  in  section  5.2.4  to  be  consistent 
with  a  simple  physical  mechanism,  that  of  vortex  pairing. 

5.2.2  Frequency  Analysis  of  Non-stationary  Signals 

In  analyzing  the  spectral  content  of  a  time  signal  such  as  that  in  figure  5.2,  the 
usual  methods  such  as  the  Fast  Fourier  Transform  can  no  longer  be  used,  as  Fourier 
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analysis  relies  on  the  basic  assumption  that  the  time  series  is  stationary.  An  alternative 
procedure  would  be  to  take  the  FFT  of  a  short  segment  of  the  data  with  the  assumption 
that  the  frequency  content  of  the  signal  is  nearly  constant  over  that  portion  of  the  data. 
This  method  has  a  serious  drawback,  though.  Taking  a  short  segment  of  the  data  implies 
a  multiplication  in  the  time  domain  with  a  suitably  chosen  windowing  function.  In  the 
frequency  domain,  this  is  then  a  convolution  of  the  spectrum  with  the  convolution  of 
the  windowing  function.  As  the  length  of  the  window  is  decreased  to  provide  a  closer 
approximation  to  a  stationary  signal,  the  effect  of  the  convolution  in  the  frequency 
domain  is  to  smear  out  the  spectral  content  and  thus  resolution  in  frequency  is  lost. 

An  alternative  method  is  provided  by  the  use  of  the  Wigner  transform.  For  the 
present  purposes,  the  Wigner  transform  will  merely  be  defined  and  then  used  to  corrob¬ 
orate  some  of  the  observations  from  the  time  series  data.  Appendix  B  will  contain  a  more 
thorough  description  of  the  transform  properties  and  the  particular  implementation  used 
in  this  study.  Additional  detailed  analyses  of  the  Wigner  transform  are  also  to  be  found 
in  the  recent  work  of  Imberger  and  Boashash  (1985)  and  Claasen  and  Mecklenbrauker 
(1980). 

The  Wigner  transform  of  a  function  g(t)  is  given  by 

w '(«,/)  =  /”»(*+  §)»'(*  -  *)«"”"’*  (5  1) 

where  g*  is  the  complex  conjugate  of  g.  In  practice  some  form  of  windowing  function 
is  also  used,  so  g(t)  would  be  replaced  by  yj(t)  =  w(t)g(t),  where  w(t)  is  the  chosen 
windowing  function.  Figure  5.3  shows  an  example  of  the  use  of  the  Wigner  transform 
on  a  typical  data  sample.  The  time  series  is  shown  on  the  left  hand  side  of  the  figure 
for  comparison.  The  standard  FFT  of  the  entire  record  is  given  at  the  bottom  of  the 
figure.  From  the  full  FFT,  one  can  see  that  the  signal  contains  most  of  its  energy 
at  the  two  dominant  frequencies,  12.5  Hz  and  25  Hz.  All  information  regarding  the 
phase,  however,  is  lost.  That  is,  we  have  no  idea  if  these  two  frequency  components 


W 

M 


were  present  throughout  the  record  or  occurred  in  some  non-stationary  manner.  The 
Wigner  transform  provides  precisely  this  information.  For  this  particular  record,  the 
Wigner  transform  was  evaluated  at  intervals  of  0.1  sec.  The  entire  spectral  content  of 
the  signal  at  each  of  these  points  in  time  is  thus  obtained.  A  Hanning  window  was 
chosen  to  minimize  the  effects  of  the  finite  record  length  of  the  signal.  This  is  especially 
important  near  the  ends  of  the  record.  The  question  regarding  the  phase  is  thus  resolved 
by  examination  of  the  time-frequency  plot  of  the  Wigner  transform.  As  can  be  seen  in 
figure  5.3,  the  energy  at  the  subharmonic  is  not  present  at  all  times.  In  fact,  substantial 
periods  are  observed  in  which  no  subharmonic  exists  at  all.  These  observations  are  of 
course  equally  well  verified  by  direct  examination  of  the  time  series  data. 


5.2.3  Comparison  with  Proposed  Resonance  Mechanisms 


The  observations  of  the  preceding  sections  have  a  direct  bearing  on  some  of  the 
current  ideas  in  the  literature  regarding  resonance  mechanisms.  As  indicated  in  the 
introduction,  a  wide  variety  of  resonance  mechanisms  between  a  two-dimensional  wave 
and  some  form  of  three-dimensional  wave  have  been  proposed  to  explain  the  onset  of 
three-dimensionality  in  boundary  layer  transition.  They  are  summarized  in  table  1.1. 
Many  of  these  mechanisms  involve  a  subharmonic  wave;  i.e.  a  wave  with  frequency  one 
half  that  of  the  fundamental  wave.  As  has  been  indicated,  the  present  data  show  the 
appearance  in  the  spectrum  of  energy  at  the  subharmonic  of  the  forcing  frequency.  It  is 
natural  to  ask  whether  this  behavior  can  be  explained  by  some  form  of  resonance.  From 
direct  observation  of  the  time  series  behavior  as  well  as  the  combined  time-frequency 
analysis  afforded  by  the  Wigner  transform,  it  seems  that  this  is  not  the  case,  since  the 
temporal  behavior  of  the  subharmonic  does  not  seem  to  correlate  well  with  that  of  the 


fundamental.  As  is  seen  in  figure  5.3,  the  subharmonic  content  of  the  signal  increases 


and  decays  without  regard  to  the  magnitude  of  the  fundamental  which  remains  relatively 


constant  throughout. 


On  the  other  hand,  it  could  be  argued  that  what  is  observed  in  figure  5.3  is  a 
non-stationary  resonance  between  the  forced  oblique  wave  and  the  randomly  occurring 
background  oscillations.  That  is,  resonance  only  occurs  when  the  naturally  occurring 
oscillations  attain  a  certain  amplitude  and  thus  is  seen  to  be  non-stationary.  This  ar¬ 
gument  implies  that  there  should  exist  a  strong  correlation  between  the  amplitude  of 
the  naturally  occurring  waves  and  the  appearance  of  non-stationary  behavior  in  the  re¬ 
sulting  time  series  data.  That  this  is  not  the  case  will  be  seen  in  section  5.3.1  when 
the  frequency  dependence  of  the  non-stationary  behavior  is  examined.  It  will  be  shown 
that  the  non-stationary  behavior  increases  with  increasing  frequency  of  the  forced  os¬ 
cillations,  and  that  there  is  little  correlation  between  the  amplitude  of  the  background 
waves  (which  are  limitted  to  a  band  of  relatively  low  frequencies)  and  the  appearance  of 
the  non-stationary  behavior  (which  is  seen  primarily  at  higher  frequencies). 

5.2.4  An  Alternative  Mechanism:  Vortex  Pairing 

An  alternative  mechanism  that  can  be  offered  as  an  explanation  for  the  observed 
behavior  is  the  simple  phenomenon  of  vortex  pairing.  If  one  were  to  move  at  the  speed 
of  the  disturbance,  the  streamlines  of  the  flow  would  appear  as  shown  in  figure  5.4(a). 
This  “cat’s  eye”  streamline  pattern  shows  clearly  that  the  disturbance  consists  of  a 
periodic  array  of  co-rotating  vortices.  In  general,  the  relatively  slow  growth  rates  of 
the  instability  waves  implies  that  the  strength  of  adjacent  vortices  will  vary  little  over 
a  few  wavelengths.  Thus,  the  induced  velocity  on  any  given  vortex  will  be  very  nearly 
zero,  as  contributions  from  those  vortices  immediately  preceding  and  those  immediately 
following  that  vortex  will  tend  to  cancel.  Any  imbalance  in  the  induced  velocity  field, 
however,  will  cause  adjacent  vortices  to  pair  together  forming  one  large  vortex.  This  is 
shown  schematically  in  figure  5.4(b).  This  type  of  vortex  pairing  is  widely  observed  in 


free  shear  flows  such  as  the  wake  and  mixing  layer,  but  has  never  been  clearly  observed 
in  the  boundary  layer. 


In  the  following  section,  it  will  be  shown  that  there  exist  conditions  for  which  the 
induced  velocity  field  due  to  neighboring  vortices  cannot  be  neglected  in  the  disturbance 
development.  The  induced  velocity  field  at  any  point  zo  is  given  by  the  Biot-Savart 
relation 


(z  -  zo)  x  Q 
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From  this  expression,  it  can  be  seen  that  the  induced  velocity  field  increases  linearly 
with  an  increase  in  the  amplitude  of  the  disturbance  (i.e.  an  increase  in  |  <5  |).  Thus,  any 
imbalance  in  the  induced  velocity  field  will  be  magnified  by  an  increase  in  the  disturbance 
amplitude.  In  addition,  it  can  be  seen  that  the  induced  velocity  field  increases  with 
a  decrease  in  the  vortex-to-vortex  spacing.  Therefore,  any  conditions  which  cause  a 
decrease  in  the  vortex  spacing  will  further  increase  the  tendency  for  vortices  to  pair. 

Several  experiments  can  be  performed  which  lend  support  for  this  proposed  mech¬ 
anism  of  vortex  roll-up.  Specifically,  the  frequency  dependence,  wave  angle  dependence, 
and  amplitude  dependence  of  the  forced  instability  waves  will  be  examined  and  will  be 
shown  to  be  consistent  with  this  mechanism.  Furthermore,  the  behavior  of  the  eigen¬ 
function  will  be  studied  to  gain  insight  into  the  vertical  nature  of  the  flow  during  the 
periods  of  non-stationary  oscillation.  Finally,  the  mechanism  of  vortex  pairing  will  be 
discussed  in  light  of  the  geometrical  structure  of  the  vortex  filaments  which  was  revealed 
in  chapter  4. 


5.3  Results  for  Oblique  Wave  Excitation 
5.3.1  Frequency  Dependence 

In  order  to  examine  the  frequency  dependence  of  the  instability  waves,  the  following 


experiment  was  conducted.  The  phased  heater  array  was  initially  programmed  to  force 
purely  two-dimensioned  waves.  The  forcing  amplitude  was  kept  constant  throughout 
the  experiment  as  were  the  freestream  flow  conditions.  Figure  5.5  shows  the  resulting 
wall-shear  time  series  as  a  function  of  the  frequency  of  the  forcing.  This  is  a  visual  cut 
through  the  stability  surface  at  a  single  Reynolds  number.  As  expected,  the  oscillations 
are  highly  damped  at  the  low  and  high  ends  of  the  forcing  spectrum,  and  in  the  center 
there  is  a  region  in  which  the  oscillations  are  most  strongly  amplified.  Figure  5.6  shows 
the  resulting  amplitude  spectra  corresponding  directly  to  the  time  series  traces  of  figure 
5.5.  The  spectra  are  seen  to  be  quite  clean;  i.e.  there  is  little  observed  response  at 
frequencies  other  than  the  forcing  frequency.  This  behavior  is  in  qualitative  agreement 
with  the  predictions  of  the  linear  theory. 

Consider  next  the  very  same  experiment,  but  with  the  phased  heater  array  now 
programmed  to  generate  a  single  oblique  normal  mode  at  an  angle  of  approximately  ten 
degrees.  The  time  series  traces  and  amplitude  spectra  are  shown  in  figures  5.7  and  5.8 
as  a  function  of  the  forcing  frequency.  A  very  striking  difference  is  noted.  At  the  lower 
forcing  frequencies,  the  observed  flow  response  is  quite  similar  to  that  observed  previously 
in  the  two-dimensional  case.  At  the  higher  end,  however,  the  response  is  quite  different. 
The  time  series  behavior  is  seen  to  be  non-stationary  with  large  amplitude  fluctuations. 
This  type  of  behavior  is  observed  to  occur  over  a  substantial  range  of  frequencies.  At 
frequencies  greater  than  40  hz,  however,  the  observed  flow  response  returns  to  the  type 
of  behavior  seen  in  the  two-dimensional  case. 

A  look  at  the  spectra  of  figure  5.8  provides  further  insight  into  this  behavior.  Ap¬ 
propriate  caution  must  be  taken,  however,  when  discussing  these  spectra  as  the  signals 
were  observed  to  be  non-stationary.  At  the  lower  frequencies,  a  very  clean  response  is 
seen  at  the  forcing  frequency.  The  response  in  this  frequency  range  is  thus  quite  similar 
to  that  of  the  two-dimensional  case.  At  the  higher  end,  however,  the  Row  response  at 


the  forcing  frequency  decreases  in  accordance  with  the  linear  theory,  but  the  level  of 
activity  throughout  the  rest  of  the  spectrum  has  greatly  increased.  This  corresponds 
to  the  observed  non-stationarity  observed  in  figure  5.7.  Beyond  a  forcing  frequency  of 
approximately  40  hz,  a  return  to  the  essentially  undisturbed  flow  conditions  is  observed. 

In  the  previous  section,  it  was  suggested  that  the  observed  non-stationary  behavior 
could  be  explained  in  terms  of  a  roll-up  of  adjacent  vortices  of  the  disturbance.  The 
frequency  behavior  observed  in  figures  5. 5-5. 8  can  now  be  discussed  in  terms  of  this 
mechanism.  Specifically,  we  must  ask  why  the  non-stationary  behavior  is  seen  to  occur 
primarily  at  the  higher  end  of  the  unstable  frequency  range.  From  equation  5.2,  it  is 
seen  that  the  induced  velocity  at  any  point  in  the  boundary  layer,  x0,  due  to  neighboring 
vortices  varies  inversely  with  the  square  of  the  distance  to  the  adjacent  vortices.  Since  it 
is  known  that  the  phase  speed  c  of  the  waves  varies  little  from  c  «  .35t7oo,  it  can  be  seen 
that  an  increase  in  the  frequency  of  the  disturbance  corresponds  to  a  decrease  in  the 
disturbance  wavelength.  Thus,  the  effect  of  increasing  the  frequency  is  to  decrease  the 
separation  distance  between  adjacent  vortices.  This  then  increases  the  induced  velocity 
quadratically  according  to  equation  5.2.  A  mismatch  in  the  induced  velocity  from  vortices 
immediately  preceding  and  those  immediately  following  a  given  vortex  will  therefore  have 
a  much  greater  effect  upon  the  tendency  for  vortices  to  roll-up.  This  mismatch  could  be 
the  result  of  random  disturbances  in  the  flow  or  perhaps  even  the  difference  in  amplitude 
due  to  the  spatial  growth  of  the  disturbance.  The  observed  increase  in  the  non-stationary 
nature  of  the  oscillations  with  increasing  frequency  is  thus  consistent  with  the  mechanism 
of  vortex  roll-up. 

5.3.2  Wave  Angle  Dependence 

A  second  parameter  that  can  be  explored  is  the  angle  9  =  tan~l(0/a)  of  an  oblique 
wave.  To  study  the  effect  of  this  parameter,  consider  the  following  experiment.  The 


forcing  amplitude  and  freestream  flow  conditions  were  again  held  constant  as  in  the 
previous  experiment.  The  forcing  frequency  in  this  case,  however,  was  fixed  at  28  Hz 
and  the  angle  6  of  the  disturbance  was  varied.  Figures  5.9  and  5.10  show  the  resulting 
time  series  and  spectra  as  a  function  of  the  wave  angle.  Two  observations  can  be  made 
from  these  figures.  First,  the  response  at  the  forcing  frequency  is  seen  to  decrease  with 
increasing  wave  angle.  This  is  in  agreement  with  the  theorem  of  Squire  (1933)  which 
stated  that  the  critical  Reynolds  number  increased  with  increasing  wave  angle,  and  thus 
the  integrated  amplification  as  recorded  in  figures  5.9  and  5.10  decreases  with  increasing 
wave  angle.  This  theorem  is  strictly  valid  only  for  temporally  growing  waves,  however. 
A  more  direct  explanation  is  afforded  by  the  calculations  of  Mack  in  which  it  is  shown 
that  the  growth  rates  for  oblique  waves  generally  decrease  with  increasing  wave  angle. 
This  is  shown  in  figure  5.11  in  which  the  integrated  amplification  ln(A/A0)  is  shown  as 
a  function  of  the  wave  angle  0. 

The  second  feature  observed  from  figures  5.9  and  5.10  is  that  the  occurrence  of  non¬ 
stationary  behavior  in  the  time  series  waveforms  is  seen  to  increase  with  increasing  angle 
just  as  was  the  case  with  increasing  frequency.  This  is  not  suprising  when  we  consider  the 
geometry  of  the  situation  as  sketched  in  figure  5.12.  It  must  first  be  recalled  from  section 
4.3  that  the  phase  speed  vector  was  found  to  be  parallel  to  the  freestream  flow  direction 
and  furthermore  was  seen  to  be  independent  of  wave  angle  9.  Thus  an  increase  in  wave 
angle  (at  constant  frequency)  leaves  A %  —  cff  unchanged  while  A,  decreases.  Therefore 
as  seen  in  figure  5.12,  the  separation  distance,  S,  between  adjacent  vortices  in  the  oblique 
wave  decreases  with  increasing  wave  angle.  The  induced  velocity  as  given  by  equation 
5.2  therefore  increases  with  increasing  wave  angle  just  as  was  the  case  with  increasing 
frequency.  This  observed  behavior  is  again  consistent  with  our  proposed  mechanism  of 
vortex  roll-up. 
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5.3.3  Amplitude  Dependence 

Returning  to  the  spectra  of  figure  5.8  (forcing  frequency  dependence)  or  figure  5.10 
(wave  angle  dependence),  an  additional  observation  can  be  made.  It  is  clearly  seen 
that  the  occurrence  of  the  non-station  ary  behavior  is  a  much  stronger  function  of  the 
forcing  frequency  or  wave  angle  than  it  is  a  function  of  the  disturbance  amplitude. 
This  is  seen  by  noting  that  at  the  most  amplified  frequency  of  approximately  24  Hx,  the 
disturbance  waveform  is  seen  to  be  quite  stationary.  At  higher  forcing  frequencies  or  wave 
angles,  however,  where  the  disturbance  amplitude  (at  the  forcing  frequency)  is  smaller,  a 
much  greater  non-station arity  hi  observed.  This  observation  is  again  in  agreement  with 
the  mechanism  of  vortex  roll-up,  for  equation  5.2  shows  that  the  induced  velocity  field 
depends  linearly  on  the  disturbance  amplitude,  whereas  the  dependence  on  frequency 
or  wave  angle  (via  vortex  separation  distance)  is  quadratic.  Thus,  the  dependence  on 
the  disturbance  amplitude  should  not  be  expected  to  be  as  great  as  the  dependence  on 
frequency  or  angle  of  the  disturbance. 

In  order  to  directly  examine  the  impact  of  disturbance  amplitude  on  the  occurrence 
of  non-stationarity,  the  following  experiment  was  conducted.  The  array  was  programmed 
to  force  a  single  oblique  wave  at  an  aagia  #  =  15*  and  at  a  frequency  of  25  Hs  The 
power  input  to  the  array  was  then  continually  increased.  Figure  5  13  shows  the  time 


series  behavior  for  this  experiment.  As  can  be  seen,  the 


to  the  forcing  m  initially 


in  accordance  with  the  linear  theory;  i.e.  a  nice  stationary  response  is  observed  At  • 
forcing  level  of  approximately  200  W,  though,  significant  portions  of  the  time  series  are 
seen  to  exhibit  the  typical  period  doubling  noa-atatioaary  behavior  that  we  have  been 
observing.  As  the  power  is  increased  further,  this  type  of  behavior  becomes  more  and 
more  prevalent  Though  no  scale  is  shown,  the  fluctuations  in  the  Anal  trace  at  times 


exceed  50  percent  of  the 


level  The  dependence  on  amplitude  shown  here  is 


not  ae  strong  as  wae  previously  seen  in  the  frequency  or  wave  angle  dependence  K*it  it 
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is  dearly  a  factor. 


5.3.4  Eigenfunction  Behavior 

k  Further  evidence  for  this  proposed  mechanism  of  vortex  roll-up  is  obtained  by  study¬ 

ing  the  vertical  nature  of  the  disturbance;  i.e.  the  eigenfunction  profile.  If  vortex  roll-up 
is  occurring,  then  the  eigenfunction  should  clearly  show  evidence  of  velocity  fluctuations 
further  from  the  plate  surface  than  would  normally  occur.  It  should  first  be  pointed  out 
that  the  occurrence  of  vortex  pairing  or  roll-up  is  almost  certainly  a  three-dimensional 
phenomenon.  Thus  the  velocity  sensor  will  see  all  three  components  of  the  disturbance 
velocity  field.  Due  to  the  high  aspect  ratio  (8:1),  though,  the  sensor  is  much  more  sen¬ 
sitive  to  velocity  fluctuations  perpendicular  to  its  long  axis,  u(y)  and  v(y),  than  it  is  to 
those  that  are  parallel,  i.e.  w(y).  Furthermore,  fluctuations  in  the  vertical  direction  are 
much  smaller  than  those  in  the  stream  direction.  The  fluctuations  seen  by  the  probe, 
then,  are  essentially  from  the  single  component,  u(y). 

In  order  to  explore  the  nature  of  the  eigenfunction,  the  following  experiment  was 
conducted  The  phased  beater  array  was  initially  programmed  to  force  an  oblique  wave 
at  15  degrees.  The  frequency  of  the  forcing  was  set  at  22  Hz.  As  can  be  seen  in  figure 
5  5,  the  response  of  oblique  waves  at  such  low  frequencies  is  very  clean.  Thus,  a  very 
stationary  oblique  wave  was  generated .  The  hot-film  velocity  sensor  then  traversed  the 
boundary  layer  in  stepe  of  0.1  mm  beginning  at  a  vertical  position  0.1  mm  from  the 
surface  of  the  plate.  Figure  5.14  shows  the  resulting  flow  response  as  a  function  of 
the  vertical  distance,  y  The  shape  of  the  eigenfunction  profile  is  quite  typical.  The 
maximum  fluctuation  level  of  umaa/lfaa  at  2%  i«  seen  to  occur  at  y  =  0  3mm  from  the 
plate 


The  experiment  was  then  repeated  with  one  small  change  in  the  forcing  conditions. 
The  forcing  amplitude  was  kept  at  the  same  level  as  before,  and  the  array  was  again  pro- 


grammed  to  output  oblique  waves  of  15  degrees.  The  only  change  was  in  the  frequency 
of  the  forcing.  The  forcing  frequency  was  set  to  30  Hz.  Reference  to  figure  5.5  shows  that 
at  this  frequency,  highly  non-stationary  behavior  is  observed.  Measurment  of  the  eigen¬ 
function  just  as  before  results  in  the  second  profile  plotted  in  figure  5.14.  The  magnitude 
of  the  fluctuations  is  seen  to  be  much  larger  than  in  the  previous  case.  This  is  not  in 
accordance  with  the  linear  theory,  as  waves  at  22  Hz  should  be  more  strongly  amplified 
than  those  at  30  Hz.  This  is  accounted  for  by  the  very  large  amplitude  fluctuations 
typically  observed  when  the  disturbance  becomes  non-stationary. 

The  important  point  to  be  observed  in  figure  5.14  is  that -in  the  neighborhood  of  the 
wall,  the  entire  character  of  the  eigenfunction  profile  has  changed.  There  is  no  longer 
a  clearly  observed  maximum.  Approximately  the  same  fluctuation  level  is  seen  from 
y  =  0.1mm  to  y  =  0.6mm.  The  large  increase  in  the  magnitude  of  the  fluctuations  occurs 
as  a  result  of  the  stretching  of  vortex  filaments  which  occurs  in  the  pairing  process.  It  is 
the  displacement  of  the  vortex  centers,  however,  that  accounts  for  the  distorted  profile 
shape  as  seen  in  figure  5.14. 

It  is  to  be  noted  that  the  maximum  fluctuation  level  in  the  case  of  vortex  pairing 
is  found  roughly  twice  as  far  from  the  wall  (y  =  0.6mm)  as  was  found  in  the  non¬ 
pairing  case  (y  =  0.3mm).  This  is  in  agreement  with  the  simplest  intuitive  notion  of 
vortex  pairing,  where  one  vortex  simply  rolls  over  another  and  thus  doubles  the  height. 
It  is  interesting  to  note  also  that  this  increase  in  activity  away  from  the  wall  has  been 
observed  in  other  experiments.  For  example,  in  the  classic  paper  of  Klebanoff,  Tidstrom, 
and  Sargent  (1962),  the  following  statement  is  made:  “The  position  of  the  maximum  in 
the  intensity  distribution  as  given  by  the  linear  theory  is  at  0.25,  but  in  the  non-linear 
range  the  position  of  the  maximum  at  a  spanwise  position  corresponding  to  a  peak  moves 
away  from  the  surface  as  breakdown  is  approached,  and  at  breakdown  has  moved  out 
to  about  0.45.”  This  is  clearly  demonstrated  in  figure  5  of  their  paper.  The  flow  in 


vW\WWvV?'v.  <<a'  \  v 

to ^ v\ ; . ; 


that  experiment  was  highly  three-dimensional  and  non-linear  having  been  disturbed  by 


a  physical  obstruction  (celophane  tape)  placed  at  regular  intervals  on  the  plate  surface. 


It  is  believed  that  the  driving  mechanism  in  that  experiment  was  precisely  the  same  as 


in  the  present  study,  namely  vortex  pairing. 


5.4  Summary  of  Three-Dimensional  Effects 


In  the  preceding  two  chapters,  the  nature  of  oblique  instability  waves  has  been  ex¬ 


plored  analytically  as  well  as  experimentally.  It  has  been  shown  that  a  great  number  of 


similarities  exist  between  oblique'  waves  and  their  two-dimensional  counterparts.  Some 


of  these  similarities  are  indeed  quite  suprising.  For  instance,  a  series  of  experiments 


involving  forced  oblique  waves  has  shown  that  several  of  the  properties  are  independent 


of  the  angle  of  an  oblique  wave.  It  was  indicated  that  the  direction  of  propagation  is 


essentially  the  same  for  all  types  of  disturbances,  two-dimensional  as  well  as  oblique. 


Specifically,  the  direction  of  propagation  is  parallel  to  the  freestream  flow  direction  in  all 


cases.  In  addition,  the  magnitude  of  the  phase  velocity  was  measured  and  was  found  to 


be  independent  of  wave  angle  as  well.  Finally,  it  was  found  that  the  region  of  unstable 


frequency-Reynolds  number  space  (i.e.  the  stability  surface)  was  also  seen  to  be  indepen¬ 


dent  of  wave  angle.  Combining  these  observations  with  our  knowledge  of  the  dominant 


forces  in  a  boundary  layer,  a  mechanism  describing  the  physical  nature  of  viscous  insta¬ 


bility  is  offered.  Namely,  it  has  been  suggested  that  viscous  instability  waves  are  caused 


by  convection  and  diffusion  of  periodic  vorticity  generated  at  the  wall.  This  explains 


why  several  of  the  observed  disturbance  properties  are  independent  of  the  wave  angle, 


as  the  mechanisms  of  diffusion  and  convectipn  occur  without  regard  to  the  orientation 


of  the  lines  of  constant  phase. 


In  addition,  it  has  been  shown  that  several  differences  exist  which  set  oblique  waves 


apart  from  two-dimensional  waves.  In  Chapter  4,  the  vorticity  fields  of  these  two  types  of 


V*N . 
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disturbances  was  examined,  and  it  was  found  that  a  very  important  fundamental  differ¬ 
ence  exists.  In  the  two-dimensional  case,  the  vorticity  field  was  purely  one-dimensional 
(i.e.  a  scalar,  u>M).  In  the  three-dimensional  case,  however,  the  vorticity  field  is  fully 
three-dimensional  with  each  component  w*,  wy,  and  w,  being  a  non- vanishing  quantity. 
It  was  then  shown  that  the  existence  of  a  three-dimensional  vorticity  field  gives  rise  to 
an  additional  source  of  vorticity  production  that  was  unavailable  in  the  two-dimensional 
case.  It  was  further  suggested  that  this  additional  source  of  vorticity  production  could 
cause  instability  to  oblique  disturbances  even  when  two-dimensional  disturbances  were 
completely  stable.  The  example  used  to  illustrate  this  possibility  was  the  as  yet  unre¬ 
solved  stability  problem  of  Couette  flow. 

Finally,  it  was  shown  that  certain  regions  of  the  parameter  space  exist  for  which 
oblique  waves  are  indeed  observed  to  be  more  unstable  than  any  two-dimensional  wave. 
This  behavior  is  believed  to  originate  from  the  additional  vorticity  production  mechanism 
outlined  above,  and  is  seen  to  manifest  itself  in  a  roll-up  of  adjacent  vortices  of  the 
disturbance  waveform.  It  was  further  shown  that  the  occurrence  of  this  vortex  pairing 
behavior  was  seen  to  increase  with  increasing  frequency  of  oscillation,  increasing  oblique 
wave  angle,  and  increasing  amplitude.  All  of  these  factors  occur  in  a  manner  consistent 
with  the  usual  Biot-Savart  law  for  the  induced  velocity  field. 


CHAPTER  6 


CONCLUSIONS 

An  experimental  study  of  the  weakly  non-linear  as  well  as  the  three-dimensional 
nature  of  boundary  layer  transition  has  been  conducted  using  a  new  technique.  The 
active  surface  heating  technique  of  Liepmann,  Brown,  and  Nosenchuck  was  employed  and 
in  this  study  was  extended  to  provide  a  means  for  controllably'and  repeatably  introducing 
three-dimensional  disturbances  into  a  laminar  boundary  layer.  A  simple  analysis  of 
the  heating  technique  was  offered  in  an  effort  to  clarify  the  mechanism  by  which  the 
fluctuations  in  the  wall  heat  flux  are  translated  into  localized  velocity  fluctuations  in 
the  vicinity  of  the  critical  layer.  It  was  shown  that  the  effective  velocity  perturbation 
introduced  by  the  heating  technique  increases  exponentially  with  the  streamwise  distance 
along  the  heater  reaching  its  maximum  value  at  the  heater  trailing  edge.  In  this  way,  the 
spatial  resolution  of  the  disturbance  is  maintained  even  when  heaters  of  large  streamwise 
extent  are  used. 

The  details  of  the  present  experimental  configuration  have  been  described.  The 
particular  heater  geometry  used  in  this  study  consisted  of  a  32-element  heater  array 
spanning  the  test  facility  at  a  single  streamwise  location.  Each  element  of  the  array  is 
independently  programmable  in  both  amplitude  and  relative  phase  making  possible  a 
wide  range  of  three-dimensional  disturbances.  The  design  and  construction  of  the  heater 
array  are  detailed  along  with  a  description  of  the  control  and  power  electronics  necessary 
to  drive  the  array. 

The  effect  of  weak  non-linearity  on  boundary  layer  transition  was  then  examined 
by  an  experimental  study  of  two-dimensional  instability  waves  of  increasing  amplitude 
It  has  been  shown  that  the  effect  of  weak  non-linearity  is  relatively  benign.  The  growth 
rates  were  seen  to  follow  the  linear  theory  up  to  perturbation  amplitudes  of  nearly  twelve 
percent  This  was  seen  by  analysing  the  local  growth  rates  as  well  as  the  integrated 


amplification  over  a  wide  range  of  frequencies.  The  only  deviation  from  the  linear  theory 
arose  in  the  form  of  a  non-linearly  generated  first  harmonic  wave  phase-locked  to  the 
fundamental,  whose  amplitude  depended  quadratically  upon  that  of  the  fundamental. 
Higher  harmonics  were  also  observed,  and  their  amplitudes  correspondingly  scaled  with 
higher  powers  of  the  fundamental  amplitude.  It  was  concluded  that  even  though  these 
non-linearly  generated  harmonic^  do  indeed  alter  the  wave  behavior  to  some  extent,  they 
are  by  themselves  not  sufficient  to  explain  the  transition  from  small  linear  oscillations  to 
the  large  amplitude,  broad-band,  three-dimensional  oscillations  characteristic  of  a  fully 
turbulent  boundary  layer. 


Furthermore,  the  relation  of  the  present  study  to  some  of  the  ideas  prevalent  in 
the  current  literature  was  discussed.  It  was  shown  that  the  results  of  the  present  study 
are  not  consistent  with  those  experimental  and  theoretical  investigations  which  report 
a  rather  low  (u/tfoo  »  1%)  threshold  amplitude  for  the  spontaneous  onset  of  three- 
dimensionality  and  the  consequent  departure  from  the  predictions  of  the  linear  theory. 
Rather,  it  is  believed  that  with  sufficient  care,  very  large  amplitude  two-dimensional 
disturbances  can  be  generated  with  little  observed  departure  from  the  relatively  slow 
growth  rates  as  predicted  by  the  linear  theory. 


The  effect  of  three-dimensionality  on  boundary  layer  transition  was  then  investigated 
through  an  analytical  and  experimental  study  of  single  oblique  instability  waves.  Single 
oblique  waves  were  chosen  (as  opposed  to  oblique  wave  pairs  or  some  other  more  elaborate 
combination  of  waves)  since  they  are  normal  modes  of  the  boundary  layer.  They  are  thus 
the  fundamental  three-dimensional  structure  in  a  transitional  boundary  layer  They 
have  the  additional  advantage  of  having  a  relatively  simple  analytic  description  which  is 
amenable  to  some  elementary  analysis. 

It  has  been  shown  that  a  great  number  of  similarities  exist  between  oblique  waves  and 
their  two-dimensional  counterparts.  Some  of  these  similarities  are  indeed  quite  suprising. 
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For  instance,  a  series  of  experiments  involving  forced  oblique  waves  has  shown  that 
several  of  the  properties  are  independent  of  the  angle  of  an  oblique  wave.  It  was  shown 
that  the  direction  of  propagation  is  essentially  the  same  for  all  types  of  disturbances, 
two-dimensional  as  well  as  oblique.  Specifically,  the  direction  of  propagation  is  parallel 
to  the  freestream  flow  direction  in  all  cases.  In  addition,  the  magnitude  of  the  phase 
velocity  was  measured  and  was  found  to  be  independent  of  wave  angle  as  well.  Finally, 
it  was  found  that  the  region  of  unstable  frequency- Reynolds  number  space  (i.e.  the 
stability  surface)  was  also  independent  of  wave  angle. 

Combining  these  observations  with  our  knowledge  of  the  dominant  forces  in  a  bound¬ 
ary  layer,  a  mechanism  describing  the  physical  nature  of  viscous  instability  has  been 
offered.  Namely,  it  has  been  suggested  that  viscous  instability  waves  are  caused  by  diffu¬ 
sion  and  convection  of  periodic  vorticity  generated  at  the  wall.  The  diffusion  is  directed 
normal  to  the  wall,  as  the  largest  gradients  occur  in  this  direction,  and  the  convection  is 
in  the  freestream  direction,  since  convection  due  to  the  disturbance  velocity  field  is  zero 
in  the  mean.  This  explains  why  several  of  the  observed  disturbance  properties  are  inde¬ 
pendent  of  the  wave  angle,,  as  the  mechanisms  of  diffusion  and  convection  occur  without 
regard  to  the  orientation  of  the  lines  of  constant  phase. 

In  addition,  it  has  been  shown  that  several  differences  exist  which  set  oblique  waves 
apart  from  two-dimensional  waves.  The  vorticity  fields  of  these  two  types  of  disturbances 
was  examined,  and  it  was  found  that  a  very  important  fundamental  difference  exists.  In 
the  two-dimensional  case,  the  vorticity  field  was  purely  one-dimensional  (i.e.  a  scalar, 
u>„).  In  the  three-dimensional  case,  however,  the  vorticity  field  is  always  fully  three- 
dimensional  with  each  component  w,,  u;v,  and  w,  being  a  non-vanishing  quantity.  It  was 
then  shown  that  the  existence  of  a  three-dimensional  vorticity  field  gave  rise  to  an  addi¬ 
tional  source  of  vorticity  production  that  was  unavailable  in  the  two-dimensional  case. 
It  was  further  suggested  that  this  additional  source  of  vorticity  production  could  cause 


instability  to  oblique  disturbances  even  when  two-dimensional  disturbances  were  com¬ 
pletely  stable.  The  example  used  to  illustrate  this  possibility  was  the  as  yet  unresolved 
stability  problem  of  Couette  flow. 

Experimentally,  it  was  shown  that  certain  conditions  exist  for  which  oblique  waves 
are  observed  to  be  more  unstable  than  any  two-dimensional  wave.  This  behavior  is  be¬ 
lieved  to  originate  from  the  additional  vorticity  production  mechanism  outlined  above, 
and  is  seen  to  manifest  itself  in  a  roll-up  of  adjacent  vortices  of  the  disturbance  wave¬ 
form.  It  was  further  shown  that  the  occurrence  of  this  vortex  pairing  behavior  was  seen 
to  increase  with  increasing  frequency  of  oscillation,  increasing  oblique  wave  angle,  and 
increasing  wave  amplitude.  All  of  these  factors  occur  in  a  manner  consistent  with  the 
usual  Biot-Savart  law  for  the  induced  velocity  field.  In  addition  the  eigenfunction  was 
measured  and  was  also  found  to  be  consistent  with  the  mechanism  of  vortex  pairing. 

Finally,  it  is  concluded  that  the  behavior  of  oblique  waves  is  much  more  compli¬ 
cated  than  was  previously  believed.  The  widely  held  belief  that  oblique  waves  are  more 
stable  and  therefore  less  dangerous  than  two-dimensional  waves  is  seen  to  be  unfounded. 
Furthermore,  it  is  concluded  that  oblique  waves  undoubtedly  play  a  major  role  in  the 
transition  to  turbulence.  It  is  to  be  noted,  however,  that  the  results  of  this  study  deal 
only  with  the  early  stages  of  transition,  and  the  resulting  high  frequency  oscillations 
characteristic  of  fully  developed  turbulent  flow  were  not  seen.  It  is  believed,  though, 
that  additional  study  on  the  further  development  of  oblique  waves  will  indeed  reveal  this 
behavior. 
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APPENDEX  A 


SIMILARITY  ANALYSIS  OP  THE  SURFACE  HEATING  TECHNIQUE 


In  order  to  obtain  an  approximate  expression  for  the  temperature  profile  above  the 
heater  element,  consider  the  situation  as  sketched  in  figure  A.l.  The  energy  equation  is 


given  as: 


V  ■  [pu(h  +  y )  -  t  •  u  +  ^  =0 


If  we  now  consider  two-dimensional,  parallel,  low  speed  flow,  we  can  make  the  following 


approximations: 


y  «  h  ,  r-u<^q  ,  u  =  (u,0,0) 


Substituting  into  A.l,  we  obtain 


dh  ,  d*T 

'“a?  “  V 


The  thermal  boundary  layer  should  be  well  within  the  linear  region  of  the  viscous  bound¬ 
ary  layer,  and  we  can  therefore  take 


r„ 

u  =  —  y 
n 


If  it  is  further  assumed  that  the  streamwise  extent  of  the  heater  is  small  compared  to 
the  distance  from  the  leading  edge,  we  may  take  rw  to  be  constant  over  the  heater. 
Substituting  A.4  into  A. 3  we  obtain 


tw  dT  d*T  k 

—  y-77  =  ,  K=  

M  0$  dy2  pcp 


If  we  now  choose  for  the  similarity  variable  r)  —  y /(a^)1/*,  and  substitute  into  equation 


A. 5,  we  obtain 


d*T  a rw  tdT 
dr\2  +  SfiK^  dtf 


(A.6) 


Taking  a  =  Z}ik/tw  and  9(r})  =  (T  -  Tw)/(Too  -  Tw),  we  arrive  at  the  simple  ordinary 


differential  equation 


0  +  „’i£  =  °  ,  9(0)  =  0  ,  »(«)  =  ! 


This  is  easily  integrated  once  to  obtain 


and  once  more,  yielding 


d$  _fi*/s 

—  =  Ce  "  /s 
dr) 


$(t j)  =C  f”  e-^^dr, 
Jo 


The  constant  in  A.9  is  evaluated  by  making  the  subsitiutions,  rf  =  (Zt)1^3  and  dr}  — 


(3t)~2/3dt  and  evaluating  at  =  oo.  Thus, 


ff(oo)  =  C  •  3~a'3  rV'e-'dt  =  C  ■  3-*/3r(i)  =  1 

Jo  ’ 


(A.  10) 


The  final  result  is  obtained  as 


${t})  =  0.776  f1  e-^^dr} 

Jo 


(A.U) 


From  this  expression  the  temperature  profile  and  the  heat  transfer  rates  above  the  heater 
are  easily  obtained. 


APPENDIX  B 


PRESENT  IMPLEMENTATION  OF  THE  WIGNER  DISTRIBUTION 

As  indicated  in  Chapter  5,  several  excellent  references  are  available  which  contain 
detailed  discussions  of  the  general  properties  of  the  Wigner  distribution.  The  present 
discussion,  however,  will  .be  concerned  primarily  with  the  particular  implementation  of 
the  Wigner  distribution  as  used  in  the  analysis  of  the  weakly  non-stationary  signals 
encountered  throughout  this  investigation.  The  Wigner  distribution  will  be  shown  to 
provide  a  viable  alternative  to  spectral  methods  such  as  the  Fourier  transform  and  will 
therefore  be  frequently  referred  to  as  the  Wigner  transform.  A  brief  introduction  to  the 
Wigner  transform  will  be  given  followed  by  several  instructive  examples  which  point  out 
several  of  the  advantages  as  well  as  the  disadvantages  of  this  method. 

The  Wigner  transform  of  a  real  function  x(t)  is  defined  as  follows: 

W{t'f)  =  /~  \)'~"vfr*  (B\) 

where  z(t)  is  an  analytic  signal  whose  real  part  is  identical  to  the  original  function  x(t), 
and  z*(t)  is  its  complex  conjugate.  The  analytic  signal  z(t)  is  defined  as 

'(,)=l(,)+;E(rV  <B21 

It  can  be  seen  to  be  formed  by  a  convolution  of  the  original  real  time  series  x(t)  with  a 
function  y(t)  given  by 

=  +  ^  (».3) 

The  Fourier  transform  of  y(t)  is  easily  shown  to  be  twice  the  unit  step  function,  i.e  : 

f  2  f>0 

Y(f)=  1  f=0  (BA) 

v  0  f<0 

Thus,  the  spectrum  *>f  z(t)  is  equal  to  twice  that  of  x(t)  for  positive  frequencies  and  is  zero 
for  all  negative  frequencies  In  the  present  implementation,  this  forms  an  operational 
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definition  of  th$  analytic  signal.  The  Fourier  transform  of  x(t)  is  calculated.  It  is 
then  multiplied  by  a  factor  of  two  for  ail  positive  frequencies  and  set  equal  to  zero  for 
all  negative  frequencies.  The  DC  component  remains  unchanged.  Taking  the  inverse 
transform  yields  the  analytic  signal,  z(t).  As  a  simple  example,  consider  the  signal, 
z(t)  =  cos(2xft).  Substituting  into  B.2  gives 

z(t)  =  e,2w^t  =  cos(2irft)  +  isin(2n  ft)  (fl.5) 

In  order  to  accommodate  signals  with  finite  record  length,  a  windowing  function 
is  also  used.  In  the  examples  to  follow  a  Hanning  window  was  chosen.  It  is  defined  as 
follows: 

,  ,  1  1  (2xt\ 

»(>)  =  j-  <B-6> 

where  Tw  is  the  width  of  the  window.  This  has  the  usual  effect  of  broadening  the 
spectral  peaks  but  is  very  useful  in  minimizing  leakage  effects  inherent  in  a  signal  of 
finite  duration.  Thus  the  transform  becomes 

W(tJ)  =  J°°"W(TMt+r-)z-(t-T-)e-'2*'rdT  (B.  7) 

In  practice  this  transform  is  evaluated  by  first  evaluating  the  kernel,  w(r)z{t  4- 
t /2)z’{t  -  t/2).  In  order  to  evaluate  this  signal  at  times  t/2,  either  interpolation  or 
over-sampling  is  used.  If  over-sampling  is  employed,  the  original  sampled  signal  x(t) 
must  be  sampled  at  a  rate  of  at  least  twice  the  Nyquist  frequency.  Standard  fast  Fourier 
transform  methods  are  then  used  to  evaluate  B.7  with  the  time  t  specified  as  a  parameter. 

In  order  to  further  explore  the  properties  of  the  Wigner  transform,  we  will  now 
take  a  look  at  several  specific  examples.  In  each  example,  the  Wigner  transform  will  be 
evaluated  analytically  according  to  B.l  and  numerically  by  the  above  mentioned  method 
using  a  sampled  version  of  the  input  signal,  x(t). 
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EXAMPLE  1 


Consider  first  the  stationary  signal, 

x(t)  =  a»n(2jr/of)  (0.8) 

The  corresponding  analytic  signal  is  easily  obtained  as 

z(t)  =  eiiwfut  (fl.9) 

Substituting  into  B.l, 

W(t,f)=  l  (fl .10) 

J -oo 

Simplifying,  this  gives 

W{t,f)  =  r  e-^-M'dr  =  6(f  -  f0)  (0.11) 

J  —  OO 

Thus  the  Wigner  transform  consists  of  a  single  peak  in  the  spectrum  at  f  ~  f0 
and  is  independent  of  time.  Figure  B.l  shows  a  plot  of  W[t,f)  for  a  sampled  version 
of  x(t)  where  fa  =  50 Hz.  The  time  series  x(t)  is  shown  on  the  left  hand  side  of  the 
plot,  and  the  standard  FFT  is  included  on  the  bottom  of  the  figure.  W(t,f)  is  evaluated 
at  intervals  of  0.1  sec  using  a  Hanning  window  of  width  0.4  sec.  The  window  is  rather 
wide,  but  temporal  resolution  is  maintained  due  to  the  nature  of  the  kernel  used  in  the 
Wigner  transform.  A  slight  ringing  is  noted  near  the  ends  of  the  record  because  of  the 
width  of  the  window  chosen,  but  otherwise  the  spectral  content  of  the  signal  is  seen  to 
be  relatively  unchanged  throughout  the  length  of  the  signal: 


EXAMPLE  2 

Consider  now  a  non-st  at  ionary  signal  of  the  form, 

x(t)  =  ain2x fit  =  ain2jr[/0  +  (/ma,  -  f0)t]t  (£.12) 
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The  frequency,  fx  -  /0  -  (/«.. 
corresponding  analytic  signal  is 


/o ft  is  a  linearly  imr*Mmn  function  *f  nm*  Th«- 


r(0  -  * 


<  0  1 1 ) 


The  Wigner  tranaform  te  calculated  exactly  as  in  the  previous  nampk  and  i*  given  «• 


"v  n  -  *•/  /i(o 


I  H  141 


That  m,  the  spectral  content  of  the  signal  at  any  time  l  cone teu  of  a  single  peak  ai 
/  =  /i  (<)  This  is  seen  graphically  in  figure  B  2  The  time  aeries  shown  imreeers 
linearly  in  frequency  from  fa  =  30 Hi  to  fm^m  =  70 Hi  As  ran  be  seen  the  standard 
PFT  shows  only  that  the  entire  eqpmJ  contains  energy  from  /  -  30 Hi  to  /  70 Hi 
All  temporal  or  phaee  information,  however,  is  loot  in  this  representation  The  Wigner 
transform,  on  the  other  hand,  shows  that  at  any  time  (.  the  signal  contains  essentially  one 
frequency  component  that  incrsaaae  linearly  with  tune  This  shows  clearly  the  advantage 
of  the  Wigner  transform,  as  both  frequency  and  phase  information  are  retained 


EXAMPLE  3 

In  the  previous  two  examples,  we  considered  signals  which  contained  only  one  fre¬ 
quency  at  any  given  time  As  a  third  sxample,  consider  a  signal  consisting  of  two 
stationary  contributions: 

x(t)  -  Mn2x/|f  *■  nn2wftt  (0  IS) 

The  analytic  signal  is  correspondingly, 

*(l)  =  (0  16) 


Substituting  into  B.l  and  simplifying,  we  obtain 
■*(*./)  =  *(/-/,)  +  *(/- A)  +  2cas2»(/l  -/,)*■«[/  - 


(A  17) 


pec  ted  The  third  term,  however  is  an  unfortunate  byproduct  of  the  non-linear  nature  of 
the  Wigner  transform  It  cowtsf  of  a  peak  located  at  the  average  of  the  two  frequencies, 
ft  and  ft  Unlike  the  other  two  peaks,  though,  it  is  modulated  in  time  by  the  factor 
coe2e(/i  /■  )C  Figure  B  3  shows  a  computed  Wigner  transform  for  a  signal  consisting 
of  two  sias  waves  at  30  Ha  aad  70  Ha  The  cross  term  at  50  Ha  appears  stationary,  but 
epoa  closer  inspection  as  shown  in  tgure  B  4,  it  w  soon  to  oscillate  in  time 

Under  certain  conditions,  thin  cross  term  can  be  avoided,  however  For  example, 
if  the  time  sense  mi  only  weakly  non  stationary  (i  s  tho  spectral  content  of  the  signal 
cheapen  slowly  with  tuns),  the  cross  term  can  be  removed  by  low-pass  filtering  W(t,  f)  in 
tuns  at  each  frequency  component  That  is,  for  sack  frequency  component  /, .  we  have  a 
tune  senas  W[t,  /,)  which  can  be  low-pass  titered  to  remove  the  cross  term  The  result 
of  this  operation  is  shown  in  tgure  B  &  in  which  the  Wigner  transform  of  figure  B  4  has 
bean  titered  with  a  cutoff  frequency  of  10  Hs  la  this  way,  noa-stationary  oscillations 
slower  than  10  Hs  are  allowed  to  remain  Higher  frequency  son-stationary  oscillations, 
however,  are  eliminated  hi  this  case,  only  the  two  stationary  components  at  30  Hi  and 
70  Hs  remain 

We  have  therefore  identited  two  types  of  time  sense  for  which  the  Wigner  transform 
provides  a  viable  means  of  spectral  analysis  The  first  type  w  a  time  series  which  contains 
energy  at  only  one  frequency  at  any  given  time  The  second  type  is  a  time  sene*  which 
may  contain  several  frequency  components,  but  whose  spectral  content  changes  slowly 
with  time  Both  of  theee  .ypes  of  signals  are  frequently  encounterred  in  transition  studies 
such  as  the  present  mvrstigation  The  observed  oscillations  are  typically  quite  band 
limittad,  and  changes  in  spectral  content  (i  s  due  to  Re  variation,  vortex  pairing,  etc  ) 
occur  rather  slowly  comparsd  to  the  frequency  of  oscillation  of  the  disturbance 


APPENDIX  C 


MICRPROCESSOR  CONTROL  AND  DATA  ACQUISITION  SYSTEM 

As  indicated  earlier,  a  small  microprocessor  system  was  used  for  control  of  var¬ 
ious  aspects  of  the  experiments  as  well  as  for  data  acquisition.  Figure  C.l  shows  a 
block  diagram  of  the  electronic  portion  of  the  experimental  apparatus.  As  can  be  seen, 
the  microprocessor  is  central  to  virtually  every  aspect  of  the  experiment.  Several  of 
the  components  shown  in  figure  C.l  have  already  been  discussed  such  as  the  sensors, 
anemometers,  and  the  phased  heater  array  control  electronics.  In  this  section,  we  will  be 
concerned  primarily  with  the  central  microprocessor  and  the  way  in  which  it  is  interfaced 
to  the  other  various  components  of  the  system. 

The  entire  system  was  designed  and  constructed  in-house,  and  it  is  the  purpose 
of  this  ssction  to  describe  only  the  principal  features  that  pertain  to  the  present  ex¬ 
periments.  The  system  is  baasd  on  an  8-bit  central  processing  unit  (CPU),  the  Intel 
8085-AH.  The  8065  has  an  6-bit  data  bus  and  a  16-bit  address  bus  enabling  the  direct 
access  of  64  kbytss  of  memory.  Though  capable  of  running  at  5  Mhs,  the  present  system 
configuration  uses  a  2  Mhs  clock  resulting  in  a  basic  clock  cycle  of  0.5 p*.  The  8085  op¬ 
erates  in  conjunction  with  an  Intel  8155  which  provides  256  bytes  of  local  “scratch  pad” 
RAM  (Random  Access  Memory),  three  I/O  (Input/Output)  ports,  and  a  timer/counter. 
In  addition,  4  kbytes  of  ROM  (Read  Only  Memory)  are  provided  on  the  CPU  board, 
the  lower  2  kbytss  of  which  contain  the  system  monitor  program. 

In  order  for  the  CPU  to  communicate  with  the  various  peripheral  devices  involved 
in  the  experiments,  a  72-bit  system  bus  was  organized  containing  the  16  address  lines,  8 
data  lines,  and  a  host  of  control  signals  established  by  the  CPU  and  the  various  pcriphals 
to  be  described.  An  8-bit  I/O  port  (Intel  8212)  was  used  to  drive  the  address  and  control 
signals  onto  the  bus,  and  a  bidirectional  bus  buffer/ driver  (8216)  was  employed  to  assert 
and/or  receive  data  on  the  data  bus.  Figure  C.2  shows  a  functional  block  diagram  of  the 
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above  mentioned  portion  of  the  system.  The  remainder  of  the  system  can  be  organized 


into  five  logical  sub-systems  as  follows: 


1.  Operator  Interface 


2.  Data  Acquisition/Storage  (via  DMA  Interface) 


3.  Data  Manipulation 


4.  Phased  Heater  Array  Interface 


5.  Stepper  Motor  Interface 


These  five  units  communicate  with  the  CPU  via  the  system  bus  and  are  schemati¬ 


cally  shown  in  relation  to  the  CPU  in  figure  C.3 


Operator  Interface 


To  provide  for  direct  access  to  the  CPU,  the  operator  interface  consisting  of  a 


keyboard/display  unit  and  an  optional  speech  synthesizer  were  made  available.  The 


keyboard  contains  24  keys  (16  hexadecimal  numbers  and  8  command  keys).  It  is  contin¬ 


ually  scanned  by  an  Intel  8279  Programmable  Keyboard/Display  Interface  which  upon 


key  closure  interrupts  the  CPU.  The  CPU  then  executes  an  appropriate  interrupt  service 


routine  which  is  contained  within  the  system  monitor.  In  this  way,  data  and/or  any  of 


the  instructions  contained  in  the  8085  instruction  set  may  be  directed  to  the  CPU. 


Data  are  displayed  by  the  CPU  by  means  of  a  6-bit  hexadecimal  display.  The 


display  elements  are  7-segment  LED’s  (Light  Emitting  Diodes)  which  are  organized  into 


a  4-digit  (hexadecimal)  address  field  (0000-FFFF)  and  a  2-digit  data  field  (00-FF)  Thus. 


the  contents  of  the  2ie  possible  addresses  may  be  examined.  An  additional  means  of 


communication  with  the  system  is  provided  by  a  speech  synthesizer.  A  Votrax  SC-01 


containing  64  possible  phonemes  is  programmed  by  the  CPU  and  interfaced  via  a  2- 


watt  audio  amplifier  to  an  8-ohm  front  panel  speaker.  This  provides  the  capability  for 


vocal  verification  of  stored  programs  and  verbal  prompts  when  data  entry  is  requested 


by  the  CPU.  In  addition,  the  vocal  output  is  quite  useful  in  alerting  the  operator  of  any 


>',V 


m 


>M.i' 


portion*  of  the  experimental  apparatus  which  might  be  malfunctioning  or  inadvertently 
left  off* line. 

Data  Acquisition/ Storage  (via  DMA  Interface) 

This  portion  of  the  system  contained  the  necessary  hardware  for  acquisition  and 
mass  storage  of  data  from  the  various  transducers  involved  in  the  experiments  The 
data  acquisition  circuitry  was  fairly  extensive  and  was  physically  housed  in  a  separate 
enclosure.  The  system  waa  built  around  the  Analog  Devices  DASI128,  a  16-channel. 
12-bit,  SO  Khs  A/D  converter.  The  16-channels  were  time  multiplexed  into  a  single 
converter,  thus  the  maximum  sampling  rate  of  50  Khs  was  achieved  if  only  one  channel 
waa  sampled.  With  all  16  channels  in  operation,  the  maximum  sampling  rate  per  channel 
is  approximately  3  Khs.  The  12  bits  of  data  together  with  the  4-bit  channel  identification 
number  are  latched  into  two  8-bit  parallel  ports  (Intel  8255)  which  are  then  transferred 
to  system  RAM  under  DMA  (Direct  Memory  Access)  control.  Upon  receipt  of  a  request 
for  a  DMA  transfer  from  the  DAS  (Data  Acquisition  System),  the  DMA  controller  issues 
a  hold  to  the  CPU.  The  CPU  then  relinquishes  control  of  the  bus  and  allows  the  DMA 
controller  to  transfer  data  from  the  DAS  to  system  RAM.  The  time  for  data  transfer  is 
therefore  much  less  than  the  time  required  for  the  CPU  to  execute  the  necessary  I  O 
instructions  for  data  transfer. 

Continuous  data  acquisition  was  achieved  by  using  a  double  buffering  scheme  Data 
were  continuously  read  into  system  RAM  which  was  divided  into  two  2K  buffers  When 
one  2K  buffer  was  filled,  data  were  then  routed  to  the  alternate  buffer  At  this  point, 
a  second  DMA  channel  would  begin  transfering  the  contents  of  the  filled  buffer  to  the 
A  ED- 6 2000  floppy  disk.  These  two  processes  of  data  acquisition  and  data  storage  were 
conducted  simultaneously  by  separate  DMA  channels  at  independently  programmable 
transfer  rates.  Transfer  to  the  floppy,  however,  was  always  faster  than  the  sampling 
thus  enabling  continuous  throughput.  Double  density,  double  sided,  eight  inch  floppy 


disks  wsrt  used  for  data  storage  and  hsd  a  maximum  capacity  of  I  Mbyte  These  were 
subsequently  stored  on  mafnetic  tape  for  further  processing  and  data  reduction 

Data  Manipulation 

This  section  of  the  system  contains  that  portion  of  the  circuitry  in  which  data  and  /or 
instructions  are  stored  and  manipulated  by  the  CPU  16  kbytes  of  RAM  are  provided 
for  temporary  storage  of  data  and  instructions,  and  8  kbytes  of  ROM  are  available  for 
permanent  storage  of  program  instructions  In  addition,  several  RAM  locations  were 
reserved  and  were  decoded  to  provide  access  to  several  peripheral  devices  on  the  system 
bus  This  technique,  known  as  memory  mapping  of  periphals,  was  used  in  conjunction 
with  I/O  mapping  in  order  to  provide  for  the  large  number  of  devices  controlled  by  the 
CPU  The  phased  beater  array,  for  instance,  was  seen  by  the  CPU  as  32  separate  address 
locations,  one  corresponding  to  each  of  the  32  heaters. 

An  arithmetic  processing  unit  (Intel  8231)  with  an  additional  4K  RAM  and  4K 
EPROM  was  also  provided  to  augment  the  data  manipulation  capabilities  of  the  CPU 
Data  were  written  into  the  APU  storage  registers  by  the  CPU  along  with  the  appropriate 
instruction  codes.  A  very  extensive  lint  of  possible  functions  (multiplication,  division, 
exponentiation,  sin,  cos,  etc.)  could  be  performed  by  the  APU  in  either  16  or  32-bit 
Axed  formats  or  a  32*bit  floating  point  data  format  As  an  example  of  the  capability  of 
the  APU.  an  PFT  algorithm  was  written  which  performed  a  1024-pomt  FFT  in  approx¬ 
imately  10  sec  The  various  capabilities  of  the  APU  were  used  for  a  wide  range  <>f  data 
manipulations,  making  the  system  as  a  whole  extremely  versatile  for  general  laboratory 
use 

Phased  Heater  Array  Interface 

A  general  overview  of  the  organisation  of  the  phased  heater  array  control  electron¬ 
ics  was  given  in  Chapter  2  It  was  shown  in  block  diagram  form  that  signals  from  a 


programmable  frequency  generator  were  fed  to  a  block  of  32  programmable  phase  and 


gain  circuits,  to  the  power  amplification  circuitry,  and  finally  to  the  heaters  The  pro¬ 


grammable  frequency  generator  consisted  of  an  8-bit  D/A  converter  (Analog  Devices 


DAC-0831)  and  a  waveform  generator  (Intersil  ICL-8038).  An  8-bit  data  word  written 


from  the  CPU  was  latched  by  the  D/A  converter  and  converted  to  an  analog  voltage 


This  voltage  level  then  determined  the  frequency  to  be  output  by  the  waveform  generator 


which  was  hardwired  to  output  a  50  percent  duty  cycle  sinewave  from  1  H*  to  50  H* 


Following  this  stage  were  the  32  phase  and  gain  control  stages  corresponding  to  each 


of  the  32  heater  elements  of  the  array.  Figure  C.4  shows  the  circuit  diagram  of  the  phase 


and  gain  control  stages.  The  circuitry  for  each  heater  was  identical  with  the  exception 


of  the  location  at  which  each  is  memory-mapped  by  the  CPU  As  indicated  earlier,  the 


CPU  recognises  the  heater  array  as  32  separate  memory  locations.  A  single  8-bit  word 


written  to  one  of  these  addresses  is  used  to  program  both  the  phase  and  gain  of  the 


heater  voltage.  The  8-bit  word  is  latched  by  an  8-bit  parallel  latch  (74ls373)  The  upper 


4-bits  then  are  level  shifted  by  a  quad  op-amp  (TL074)  making  the  levels  compatible 


with  a  quad  CMOS  analog  switch  (4046).  This  switch  is  used  to  selectively  bypass  any 


of  the  4  resistor  values  (10K,20K,40K,or  80K)  which  form  part  of  a  constant  amplitude 


phase-shift  stage.  This  provides  16  discrete  phase  shifts.  The  output  of  this  stage  is  then 


directed  to  the  next  neighboring  phase-shift  stage  (i.e  the  next  heater)  and  to  the  gain 


stage  as  well. 


The  lower  4-bits  follow  a  similar  path  through  a  level  shift  to  an  additional  analog 


switch.  In  this  stage,  the  switch  selectively  connects  the  output  of  the  previous  phase- 


shift  stage  to  any  or  all  of  the  4  resistor  values  (20K,40K,80K,160K)  which  form  part  of 


a  summing  amplifier.  If,  for  example,  the  signal  was  connected  to  all  of  the  resistors. 


the  equivalent  impedence  of  this  parallel  configuration  would  be  10K  (the  same  value 


as  the  feedback  resistor),  and  this  stage  would  therefore  function  simply  as  a  unity  gain 


amplifier.  If  on  the  other  hand,  the  signal  from  the  phase-shift  stage  was  connected  only 
to  the  160K  resistor,  then  the  output  voltage  would  be  1/16  of  its  original  value  tn  this 
way,  16  discrete  voltage  levels  are  possible.  The  final  stage  provides  a  trimming  amplifier 
which  is  used  to  compensate  for  non-uniformities  in  the  heaters. 

Stepper  Motor  Interface 

The  final  section  of  the  digital  hardware  which  is  important  in  the  present  study 
is  the  connection  to  the  mechanical  world,  namely  the  interface  to  the  boundary  layer 
traversing  mechanism  via  a  feedback  controlled  stepper  motor.  A  4-pole  stepper  motor 
with  a  resolution  of  200  steps/revolution  was  used.  The  motor  was  driven  by  a  set  of 
4  appropriately  phased  pulse  trains  which  were  generated  by  a  programmable  interval 
timer  (Intel  8254).  The  pulses  were  converted  to  a  compatible  analog  signal  by  4  D/A 
converters  (MC-3410).  An  optical  encoder  (Datametrics  K15  series)  was  mounted  on  the 
shaft  of  the  stepper  motor  and  provided  feedback  information  on  the  exact  position  of  the 
traverse.  The  output  of  the  encoder  was  connected  to  the  CPU  by  way  of  2  of  4  available 
A/D  converters  (Analog  Devices  AD570).  Both  the  A/D  and  the  D/A  converters  were 
interfaced  to  the  CPU  through  several  Programmable  Peripheral  Interfaces  (Intel  8255) 

The  stepper  motor  was  then  connected  to  a  traversing  gear  through  a  flexible  shaft 
coupler.  The  traverse  converted  one  revolution  of  the  shaft  into  a  10  mm  displacement 
normal  to  the  plate  surface.  During  boundary  layer  profile  measurements,  the  increment 
in  the  vertical  direction  wae  0.1  mm  or  one  tenth  of  a  revolution.  The  hot-film  velocity 
probe  was  connected  to  the  traverse  by  the  supporting  structure  as  discussed  in  Chapter 
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Figure  3.7  Amplitude  of  First  Harmonic  vs.  Forcing  Level 
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Figure  5.3  Wigner  Transform  of  a  Non- Stationary  Time-Senes 


FIfur*5.4  Alterred  Streaallne  Pattern  Due  to  Vortex  Pairing 
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Figure  B.2  Wigner  Transform  of  a  Non-Stationary  Signal 
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